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Automated Visual Inspection 
for Printed Circuit Boards 


@ Moritoshi Ando @ Masato Nakashima @ Takefumi Inagaki 


(Manuscript received March 4, 1988) 


Fujitsu has been developing automatic visual inspection systems for printed circuit boards 
(PCBs) since 1979. This paper describes newly developed inspection technologies for inter- 
mediate products of PCBs such as photomask patterns, inner layer copper patterns, plated 
through holes of multilayer boards, and mounted components. New techniques were 
developed in such areas as sensing optics and inspection algorithms. 

The systems developed have been installed at Fujitsu’s plants and have contributed to the 


high reliability of PCBs. 


1. Introduction 

Production lines for manufacturing printed 
circuit boards (PCBs) have mostly been 
automated except for visual inspection. This 
inspection has conventionally been done with 
the naked eyes using a magnifier, and it has been 
very difficult to automate. 

Human inspectors have difficulties in keep- 
ing stable and reliable inspection. The miniaturi- 
zation of printed wiring board (PWB) patterns 
used in computers has made it difficult for 
humans to do an efficient job of inspecting these 
products. It is thus necessary to develop 
automatic inspection techniques to insure 
product reliability. 

Image data processing has been undergoing 
dramatical developments over last ten years, 
and some automated inspection systems have 
been developed. But there are still some 
remaining difficulties in image data processing 
when used in practical systems. The reason 
is that it is difficult for image data processing 
to be accurate when evaluating a wide variety 
of criteria such as color, surface roughness, 
and light reflectance of the inspection targets. 
Therefore, automization of inspection systems 
depends heavily on image capture techniques 
that are able to accomplish this. 
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This paper describes newly developed 
sensing methods and inspection algorithms, 
that make automatic systems possible. 


2. Overviews of development 
2.1 Roles of inspection for PCB production 

Figure 1 shows the process of PCB produc- 
tion. Visual inspection is performed at each step. 
In the production line, photomasks may have 
several defects, and inner layer sheets may 
also have some defects. These intermediate 
products must be inspected and repaired before 
proceeding to the next process. Here, the PCBs 
are finished products where components are 
mounted on PWBs. 

Figure 2 shows an evaluation of the inspec- 
tion cost at each stage of PCB production. 
As with any manufacturing process, the cost 
of defects increases as items get closer to the 
finished product. In the case of PCB production, 
it is also remarkably cost effective to inspect 
at early stages and remove the defective pieces 
as soon as possible. 


2.2 History of development 

Fujitsu developed an automatic photomask 
inspection system in 1979, the first of its 
kind. This system used laser scanning optics 
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Fig. 2—Economical evaluation of inspection. 


and a specialized algorithm to detect non- 
permissible defects only. 

Figure 3 shows the history of development 
as it progressed from inspecting photomasks, 
until it reached the final stage of inspecting 


mounted components on the PCBs. The inspec- 
tion technology becomes more complex as the 
objects approach the finished products. For 
example, photomask inspection handles two- 
dimensional binary image data, and inspection 
for mounted components uses three-dimensional 
data. 


3. Key technologies 
3.1 Sensing and optical processing 

Image capture of PCBs should be performed 
at a speed of more than 40 MHz, and at a 
contrast greater than three. 

Table 1 lists the sensing types used. Light 
sources and detectors are chosen according 
to the object. Special optical processing, such 
as filtering, masking, and combination optics, 
have been developed for high contrast image 
capture. 
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Table 1. Sensing types 


Objects Sensing data types Light source Detector, process 
Photomask Spot light transmittance He-Ne laser Photomultiplier 
Inner layer Reflectance (3-D SCAN) He-Ne laser Photomultiplier + spatialfilter 
Plated through hole Leakage-light transmittance Halogen lamp CCD linesensor + roller mask 
Mound components Structured light reflectance He-Ne laser CCD areasensor 
Diffused-light reflectance Halogen lamp | CCD areasensor 


1978 1980 1982 1984 1986 1988 


Photomask 12/1979A 
Inner layer 2/1983 A 


Plated through 12/1984 A 
hole 

Mounted —I/1985 A 
components 


Fig. 3—History of development. 


3.2 Inspection algorithms 

Inspection algorithms have two require- 
ments, which are necessary to make the systems 
practical. 

1) Non-reference inspection 
2) Non-permissible defect detection 

The former is for easy inspection prepara- 
tion. Algorithms for non-reference is based 
on pattern rule judgement and pattern feature 
extraction. 

The latter is necessary for easy repairs of 
defects. Algorithms for non-permissible defect 
detection is based on the combination of pattern 
width measurement, which detects pattern 
defects, and feature extraction, which detects 
random defects, and together they detect non- 
permissible defects selectively. 


4. Photomask pattern inspection 

If a photomask has any defects, they are 
transcribed onto the PWBs during the photo- 
etching process. Production of the photomask 
free of defects is an important requirement for 
improving the quality and production efficiency 
of PWBs. 
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Inspections and repairs are necessary steps 
in the production process, because each photo- 
mask may have several defects. The pattern is 
generated with a computer-aided-design system 
(CAD), and then developed on photographic 
film. Scratches on the film or dust in the 
solution can cause most of the pattern defects 
with sizes ranging from ten to several hundred 
microns. Sometimes there are several hundreds 
of defects in a single mask. 

These patterns have always been inspected 
and repaired by human inspectors. However, 
as photomask patterns have become finer and 
more complex, it has become difficult for 
inspectors to inspect them. Therefore, an 
automatic inspection method is required. 


4.1 Requirements 

An automatic photomask pattern inspection 
system must be capable of sensing the pattern 
with a high S/N ratio, and discriminating 
between permissible and non-permissible 
defects. Details of these requirements are 
described in the following section. 

4.1.1 Requirements for sensing 

Only patterns that cause defective copper 
patterns need be detected. For example, dust 
on a black pattern does not cause any defects 
because the photoetching process uses _trans- 
mitted light. However, dust on a transparent 
area causes excessive defects on the PWBs. 
Scratches do not result in defects when they are 
on the transparent part. However, scratches 
in a black pattern cause open areas on the PWBs. 
Therefore, the sensing optics should be able to 
detect dust located on transparent areas and 
a scratches located on black patterns. 
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Fig. 4—Non-permissible defects. 


4.1.2 Requirements for inspection algo- 

rithms 

The inspection algorithms should detect 
only the non-permissible defects. Non-permis- 
sible defects are defined as those that damage 
the performance or quality of the PWBs. Permis- 
sible defects are ones that do not effect the PWB 
quality. In other words, permissible defects 
are either ones that are located in free space 
areas where defects do not damage the PWB 
patterns or those that are so small in size that 
they disappear after the photoetching process 
is completed. 

Figure 4 shows examples of non-permissible 
defects. Figure 4-a) shows a mouse bite with 
a depth of 100 um, and Fig. 4-b) shows a 
protrusion with a height of 80 um. 

Figure 5 shows the critical areas for judging 
whether a defect is permissible. The mask 
surfaces are divided into four inspection areas; 
lead area, pad area, inside pad area and free 
area. The critical areas are those within the 
pattern and up to 100 um away. Table 2 
indicates defect sizes and types in the critical 
area. The sizes of non-permissible defects are 
defined according to their locations and types. 
Smaller defects have to be defected in the 
lead area, rather than other areas. Their permis- 
sible sizes also depend on their types. For 
example, shorts and breaks are always non- 
permissible, as are protrusions and mouse bites 
more than 30 wm wide. In the free area and 


inside pad area, however, large defects are 
permissible. 

Thus, the sizes of non-permissible defects 
are determined according to the pattern area 
and the types of defects. 


4.2 Sensing method 

A laser beam scanning method is used to 
detect the pattern. Transmitted light intensity 
is sensed by a photodetector. The configuration 
of the system is shown in_ Fig. 6. A He-Ne laser 
beam with a 633 nm wavelength is focused 
by a scanning lens and scans over a mask by 
means of a rotating polygon mirror. The laser 
beam diameter on the mask is 10 um at half 
the maximum intensity of the Gaussian distribu- 
tion. The laser beam scans 50mm in the 
direction of the x-axis. The effective read length 
is 41mm. The mask travels in the direction 
of the y-axis on an X-Y table in 10 um steps 
and is continuously scanned. 

A photodetector receives the beam coming 
through the mask to obtain pattern signals. The 
laser beam also scans over a grating mask in 
order to detect the location of the scanning 
laser beam. These signals are used to control 
the scanning speed. 

The pattern signals are sampled at 10 um 
per bit (size of picture element) with the aid 
of the grating signals. 

The laser beam scanning system provides 
the following advantages: 
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Table 2. Defect sizes and types in critical area 


Pattern width Random 
avaed-tnes Mouse-bite Shorts Breaks Black spots 
: yp Protrusion Pinholes 
Locations a Bccess 
Lead area 230 um =10 um 230 um 
Critical area Pad area = 80 um 210 um 230 um 
Inside pad area 220 um 280 um 
Free area = = 80 um 
7 Pattern signals 
Ss Analog & digital 
Permissible * Lea, Ta giants 
defects Photodetector 
100 «wm 
[ Grating 
H signals 
| Lead area 
| : Critical 
f Pad area areas 
| 
Free Inside pad area 
! 


Fig. 5—Critical areas and defects. 


1) Spot illumination enables pattern shapes 
to be read with ease, that is, the S/N ratio is 
high 
Using transmitted light for the sensing 

system has the advantage that dust on black 

patterns and scratches on transparent film do 
not affect the pattern shapes sensed. 

2) Long focus enables pattern to be read clearly 
A focus depth of +300 yum is obtained. 

This gives sufficient depth for the inspection 

system. 


4.3 Inspection algorithms 

4.3.1 Algorithm to detect non-permissible 

defects 

The detection logic consists of pattern 
width measurement and the Sensing Matrix 
Region Tracing (SMART) method. The pattern 
width measurement detects the widths of the 
patterns and SMART detects random defects 
such as pinholes and black dots. It is necessary 
to combine these methods to detect most 
of the non-permissible defects. 
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4.3.2 Pattern width measurement 

Automatic measurement of pattern width 
without artwork data for comparison requires 
the directions and the areas to be measured. 

An algorithm using the parallel edge detec- 
tion method has been devised to obtain this 
information directly from the patterns to 
be inspected. 

Figure 7 shows how pattern width measure- 
ment is done. 

1) The edge signals in four directions (0°, 45°, 
90°, and 135°) are extracted from the 
pattern signals. 

2) After detecting some sets of parallel edges 
from the patterns in the recognition 
windows which have 64 x 64 bit sensors, 
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Fig. 7—Pattern width measurement procedure. 


measure the pattern width perpendicular 

to the edges. 

3) Compare the measured results with the 
standard ranges of pattern width, for 
exmple, W, +a and W, +8. Where W, and 
W, represent standard widths and a and 6 
the tolerant ranges. A pattern width defect 
is detected when the measured result is 
outside these ranges. 

This method can find the directions and 
areas to be measured directly from the inspected 
patterns. Consequently, when two or three 
sets of standard widths and their tolerant ranges 
are input, the inspection system automatically 
sets the standard ranges for pattern width. 
Figure 8 shows the areas inspected by the 
pattern width measurement method. 


Fig. 8—Areas inspected by pattern width measurement. 


4.3.3, SMART detection 

The SMART (Sensing Matrix Region 
Tracing) method was developed to detect 
defects within a certain distance inside or 
outside pattern edges. 

Figure 9 shows the inspection areas and 
matrices of detection bits of SMART. As 
specified number of detection bits are lined up, 
beginning from the pattern edge inside pattern 
area, or outside of the pattern within the critical 
areas. This SMART sensor is set in eight direc- 
tions, each separated by 45 degrees. 

The SMART method detects random defects 
along normal pattern edges and_ inspects 
clearances between neighboring patterns. The 
pattern detection window is 16 bits square. 
Figure 9-b) shows four groups of the sensors 
employed. 

This method can also be applied to detect 
the space between adjacent patterns. 


4.4 System configuration 

Figure 10 shows the automatic mask inspec- 
tion system for the PWB photomask. The 
system consists of an inspection optical unit, 
a detection logic unit, an inspection result 
output unit, and a system control unit. 
Figure 11 shows a photograph of the system. 
The main specifications are listed in Table 3. 
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Fig. 9—Areas inspected by SMART detection and SMART matrices. 
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Fig. 10—Block diagram of the inspection system. 


4.5 Inspection characteristics 

Figure 12 shows the detection rate for non- 
permissible defects comparing conventional 
inspection and automatic inspection. It shows 
that SMART detects 74 percent of the non- 
permissible defects and pattern width measure- 
ment detects 57 percent of them. 

By using the pattern width measurement 
method and the SMART method together, 
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Table 3. Main specification of the photo mask 
inspection system 


ee 
Item | Specifications 


Resolution 10 um 


Inspection time 3 min/300 mm square 


Defect map 
Out 
a (With 2.5 mm square marks) 
Size 15% 1.5% 12mm(@ % 6 xh) 


Fig. 11—Automatic mask pattern inspection system, 


all the critical areas of the mask can be 
inspected. Pattern width measurement used 
alone detected 100 percent of the pattern 
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Fig. 12—Detection efficiency for non-permissible detects of the conventional 
inspection and automatic inspection. 


width defects, and SMART used alone detected 
100 percent of the random defects. Combining 
these methods allowed all the non-permissible 
defects to be detected. 

Figure 12. shows a comparison between 
detection efficiencies of non-permissible defects 
for visual inspection and the new inspection 
system. The detection rate of visual inspection 
is about 70 percent for negative masks and 
55 percent for positive masks. These results 
indicate that visual inspection alone is not 
sufficient; two or three repeated inspections 
are required to obtain a detection rate of at 
least 90 percent. Figure 12 shows these results. 
Visual inspection detects many permissible 
defects; this is shown by the large circles. 

The automatic inspection system is capable 
of detecting 100 percent of non-permissible 
defects for both negative and positive film 
masks. 


5. PWB inner layer pattern inspection 

High reliability is required for the inner 
layers because they are laminated to make 
multilayer boards. Defects on the inner layer 
patterns must be detected and repaired before 
the lamination process. 


The inner layers are usually inspected 
visually and with an electrical continuity tester. 
The tester can only detect shorts and breaks 
in the conductive pattern, but cannot detect 
mouse bites and protrusions which affect 
reliability. Conventionally, visual inspection 
with an optical magnifier is combined with the 
tester to detect all of the defects. 

However, the widths of the conductive 
patterns in newly designed boards are often 
less than 100 um, and a minimum resolution 
of less than 10 wm is needed. It is very difficult 
for humans to inspect such fine patterns effec- 
tively. 

To solve these problems, a new automatic 
inspection system has been developed that 
can inspect inner layer conductive patterns 
using an optical method. 


5.1 Requirements 
The following three requirements must 
be met before the inspection system can become 
practical. 
1) Detection of both the top and bottom width 
The conductive patterns are formed on 
epoxy or polymide resin substrate with the 
etching process. The pattern has a 100 wm width 
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and 35 um thickness. The cross section is 
a trapezoid. To assure conductivity of the 
pattern, mouse bites, poor spacing and protru- 
sions must be detected in addition to shorts 
and breaks. The inspection of both the top 
and bottom width of the conductive pattern 
is required. 

Mouse bites usually have a narrower pattern 
width, but they sometimes have the correct 
bottom width and too thin a thickness. To 
detect this type of mouse bite, as shown in 
Fig. 13, it is necessary to measure the top width 
of the conductor. 

Poor spacing sometimes has the correct top 
width, but a defective bottom width as shown 
in Fig. 14. In this case, it is necessary to measure 
the bottom width of the conductor. 
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A CCD line sensor is used as the sensing 
device in conventional systems. The top width 
of the conductive pattern can be detected, 
but the bottom width of the pattern cannot 
because the side wall reflects poorly. 

2) Detection of all defects 

Defects are divided into two types. The first 
type of defect is found in wide areas such 
as incorrect width, and the second type is 
random defects such as mouse bites, protrusions, 
and opens. The algorithm must be able to 
extract all of these defect types. It must also be 
able to filter out permissible defects such as 
small isolated patterns or dust patterns that 
do not affect continuity. This is important 
because if these permissible defects were 
detected, the yield would decrease. 
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Fig. 15—3D SCAN detection. 


3) Easy indication of defect positions 

Defects such as excess copper can be 
repaired by a skilled worker. Easy-to-use output 
is required so that the defects can be easily 
found and repaired. 


5.2 Sensing method 

5.2.1 Conditions of pattern detection 

The following conditions are required. 

1) Top and bottom width of copper pattern 
should be detected. 

2) Accuracy should be better than 10 wm. 

3) Image contrast should be larger than three. 

4) Depth of focus should be larger than 
+300 um. 

5.2.2 3D-SCAN detection 

3D-SCAN detection method was developed 
to meet the above requirements. A laser scanner 
is used and both the top and the bottom width 
of the copper line can be detected. The configu- 
ration is shown in Fig. 15. The He-Ne laser beam 
is deflected with a rotating polygon mirror 
and is scanned on the inner layer pattern. 
Deflected light from the inner layer is retro- 
reflected and imaged on a plane by a lens. 
A spatial filter with a pinhole of diameter 
100 um in the center is set on an imaging plane. 
The transmitted light through the pinhole is 
detected with a photomultiplier tube. 

When the beam strikes position A, it diffuses 
into substrate material before it is reflected. 
The image of beam A’ on the image plane 
increases in size and has a flat distribution. 


10 


metres hold-for bottom width 


__ 


Fig. 16—Copper pattern signal. 


When the diameter of the filter is set to the same 
size as the scanning beam, the filter blocks 
the reflected light beam from the substrate. 
However, when the laser beam strikes the 
top of conductor C or side wall B, it does 
not diffuse and the image sizes of laser beam 
B’ and C’ are the same as the scanning laser 
beam. When the filter size is larger than the 
scanning laser beam, most of the light energy 
passes through the filter. This method enables 
distinguishing the substrate from the top of 
the copper line and side wall line with a suf- 
ficient S/N ratio. 

The result is shown in Fig. 16. Large sig- 
nal ‘“‘a” 
flected light signal from the top of the cop- 
per line. Small signals ‘“‘b” at the sides of main 


at the center corresponds to the re- 


signal ‘“‘a” correspond to the signals from the 
side walls of the copper pattern. Here, two 
threshold levels are used to measure the top 
and bottom width of the conductive pattern 
simultaneously. 

Defects detected using this method are 
shown in Fig.17. Both the top width area 
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Fig. 17—Examples of defects. 
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Fig. 18—Stage height deviations vs. measured copper 
pattern width. 


(white) and bottom width area (gray) can be 
detected. 
Figure 18 shows the focal depth obtained 
which was +300 wm with an accuracy of +5 wm. 
This system can obtain a deeper focal depth 
than conventional CCD camera systems. 


5.3 Inspection algorithms 

The algorithms are a combination of the 
width measurement and feature extraction 
methods. They can detect both pattern width 
defects and random defects. 
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Fig. 19—Pad width measurement method. 


5.3.1 Width measurement 

The width of the conductor lines are 
measured to detect such items as defective 
pattern width and poor pattern spacing. 

PWB patterns consist of both leads of 
100 um width and pad patterns of 800 um 
size. Two algorithms have been developed 
for both pad and lead inspection. 

5.3.2 SMART detection 

The SMART (Sensing MAtrix Region 
Tracing) can inspect whether the pattern width 
satisfies the tolerance. The working principle 
is the same as that developed for photomask 
inspection as shown in Fig.9. This SMART 
sensor is set in eight directions, each separated 
by 45 degrees. All pattern gap defects can 
be also detected when the same sensors are 
applied to the pattern space. 

5.3.3. Pad measurement 

This algorithm is used to measure the 
diameter of the pad. The pads are usually 
rectangular or circular and are sometimes 
hexagonal. When rectangular shape is inspected, 
the spacing between two pairs of parallel edges 
are detected and measured. When a hexagonal 
shape is inspected, the spacing between three 
pairs of parallel lines are measured. An example 
of the sensor for inspecting a hexagonal shape is 
shown in Fig. 19. A to H indicate line measure- 
ment sensors. 
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Fig. 20—Local defect detection algorithm. 
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Fig. 21—Automatic PWB pattern inspection system. 
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The inspection procedure is as follows. 
1) Pad center detection 

Pad center is detected when the results of 
sensor A to C and B to D are equal. 
2) Parallel edge measurement 

The results of sensor E and G, and F and H 
are added. These pairs are parallel to A and B. 
The same value as ones at the center would 
be obtained if the hexagonal shape is normal. 

3) Defect detection 

Measured results are compared to the 
tolerance value. It is determined whether the 
results are within the tolerance range. Measure- 
ments are done only at the center area. The 
four sensor pair measure the whole area and 
the results are then compared to the tolerance 
values. 

In the example of Fig. 19, the pad pattern 
has a mouse bite. When added values F and H 
become less than the standard values, the pad 
is determined as defective. 


The measurements of a rectangular and 
a circle pad can be done using these sensors. 

5.3.4 Local defect detection 

The local defect detection algorithm can 
detect random defects such as breaks, protru- 
sions and extraneous coppers. It detects abrupt 
changes in edge orientation which occur at 
the defect. Figure 20 shows a diagram of the 
algorithm. Sensor A is a copper pattern detec- 
tion sensor and sensor B is a substrate detection 
sensor. When the each sensor detects the cor- 
responding pattern, the pattern is determined 
as defective. For example, Fig. 20-a) shows 
an example of a break detection and Fig. 20-b) 
shows an example of a protrusion pattern. 
Breaks and protrusions in all direction can be 
detected by using eight sensors each 45 degrees 
apart. These sensors inspect both the top and 
the bottom width pattern. They can detect 
not only breaks and protrusions of the bottom 
width pattern, but can also detect mouse bites 
in the top width. 


5.4 System configuration 

Figure 21 shows the automatic inspection 
system, and Fig. 22 shows its configuration. 

The system is made up of the following 
components. 
1) System controller 

The system controller controls X-Y stage 
movements, memorizes the positions of defects 
and controls the threshold values that dis- 
tinguish between defective and normal patterns. 
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Inspection algorithms 
Binary pattern 


He-Ne laser 


Defect position output 
Pattern detection 
unit 


Rotating mirror 


Stage controller 


Fig. 22—Schematic diagram of inspection system. 


2) Pattern detector 

The pattern detector consists of a laser 
scanner and an X-Y table as shown in Fig. 23. 
The laser beam is guided to a twelve faceted 
mirror rotating at 7 200rpm. The resolution 
is set to 10 um, and the scanning speed is 
100 m/s. Reflected light is detected with the 
3D-SCAN detector, and both the top and 
bottom widths are detected. The PWB is firmly 
attached to the X-Y table with a vacuum. 
A PWB 300 x 360 mm in size can be inspected 
in three minutes. 
3) Defect detector 

The detector to detect defects consists of 
two width measurement algorithms, one for 
line patterns and the other for pads, and two 
feature extraction algorithms, one fora SMART 
sensor and the other for a protrusion sensor. 
In this system, all inspection algorithms are 
implemented in hardware for high speed. 
4) Detected results display 

The system has two output devices, a CRT 
display and a hole punch. The CRT displays 
defect shapes. The defective positions are 
punched on a paper sheet the same size as the 
inspected PWB. This is shown in Fig. 24. The 
holes are 2.56 mm in diameter. The paper is 
laid over the PWB, and the defects can be seen 
through the punched holes. This allows the 
defects to be easily found. The defects are 
examined to determine whether they can be 
retouched by a human worker. The repairable 
defects are retouched through the punched 
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Fig. 23—Optical apparatus. 


Punched holes 


Alignement holes 


Fig. 24—Defect location reporting method by punched 
holes. 


holes. The inspection results for both sides of 
a PWB are punched out on two consecutive 
pages. The pages are then wrapped around the 
PWB, and also protect the PWB from scraches. 


5.5 Inspection characteristics 

The main characteristics of the system 
are listed in Table 4. This system was imple- 
mented at a plant and evaluated. Figure 25 
shows the results. PWBs were inspected by both 
the automatic inspection system and the con- 
ventional visual inspection. It shows that the 
system can detect all types of defects. Visual 
inspection yielded relatively good results when 
detecting breaks, shorts, and extraneous cop- 
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Table 4. Main characteristics of the inner layer 
inspection system 


Items Specifications 


Resolution | 10 um 


Inspection time 3 min/300 mm square 


Criteria values 
Output 


Independently selectable 


Punched hole sheet (dia 2.56 mm) 


Excessively wide copper 


Mouse bite -—7~ | ie Too narrow 
f Automatic 7 \\ 
/ inspection \ 
/ ' \ 
: \ 
. \ 
y 


M = xy Poor spacing 


Human 
inspection | 


Extraneous copper 


Fig. 25—Comparision of performance between 
automatic and visual inspection. 


pers, however, missed patterns that were ex- 
cessively wide or too narrow and also missed 
poorly spaced pattern. Human _ inspectors 
also missed protrusions and small mouse bites. 

The system shortened the inspection time 
down to 4 min from what took human inspec- 
tors 20 min to 40 min, and could detect not 
only bottom width defective patterns but 
also top width defective patterns. 

As for inspection speed and resolution, 
further development has been continued by 
increasing image processing speed and using high 
accuracy optics. Up to now, inspection speed 
has been increased ten fold with 5 um resolution 
on the newly developed system. 


6. Through hole inspection 

PWB technology is getting more complex 
as LSI technology developes. As a result, the 
number of inner layers of the PWB are increasing 
and the diameter of the through hole, smaller. 
The PWB used for the FACOM M-380 for 
example, has through holes with diameters 


0.2 mm Crack Void 


| 


iu 


! 


I 


Plated solder Substrate 


Fig. 26—Cross section of PWB. 


as small as 0.2 mm, and aspect ratios as large 
as 15. Plated through holes connect electrical 
circuits of inner layers vertically through the 
substrate with electrodepositioning copper on 
the hole wall. 

Through holes sometimes have defects 
such as voids and cracks as shown in Fig. 26. 
Electrical continuity testers and visual inspection 
are the conventional means of inspecting 
through holes. The electrical continuity tester 
can detect plating voids that break the continui- 
ty of the circuit, but they cannot detect cracks. 
Visual inspection is effective in detecting cracks 
only in holes with low aspect ratio PWBs. It 
is not efficient for inspecting high aspect ratio 
through holes. The automatic inspection system 
for the high aspect ratio through hole is very 
much required. 


6.1 Requirements 

The requirements for inspecting plated 
through holes are high speed and non-reference 
inspection. The system must be capable of 
inspecting a 300 mm square PWB within a few 
minutes. Non-reference inspection is also 
required because large scale memory is required 
to store the hole position data. 


6.2 Sensing method 

6.2.1. Principle of sensing method 

A new method called leakage light detection 
has been developed using the optical methods as 
shown in Fig. 27. 
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Defective through hole Normal through hole 


Leakage light 0.2 Ix 


Illuminating light 300 kIx 


Fig. 27—Leakage light detection. 


Defective through holes Normal 


Fig. 28—Leakage light from defective through hole. 


The method uses following four facts. 
1) Plated through holes must be made of metal 
cylinders. 
2) The substrate is exposed to the inside of the 
plated through hole. 
3) The substrate diffuses and transmits some 
light. 
4) The light can pass through the defects. 
The detection process is explained as 
follows. 
One side of the hole is covered with a mask. 
When the hole is illuminated with intense light, 
the light diffuses into the substrate. When it 
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Discrimination algorithm 


Substrate transmitted light 


Warpage 


Noise light 


Leakage light 0.2 Ix 


Illumination 300 klx 


Complete masking 


Fig. 29—Technical requirements. 


reaches a defect, the light is transmitted through 
the defect, and leakes into the through hole. 
When the PWB is observed from the other side, 
only defective through holes can be seen, while 
normal through holes remain dark. Then this is 
the ‘Leakage light detection” method. 
Figure 28 shows the experimental results. The 
three bright holes at the left indicate defective 
through holes, and the two dark holes at the 
right show the normal through holes. 

6.2.2 Leakage light detection 

There are three technical requirements to 
make the method practical as shown in Fig. 29. 
1) Detection of weak light 

Leakage light intensity from detects is low. 
For example, when PWBs are illuminated with 
300 klx, the leakage light is only 0.2 lx. The 
system has to be able to detect this weak light. 
2) Complete masking 

The through hole being inspected must be 
completely masked. This is so the illuminating 
light will not leak into the through holes. If it 
does leak, the S/N ratio becomes worse. There- 
fore the masking has to be complete in spite of 
PWB warpage as large as | mm. 
3) Discrimination between defective leakage 

light and substrate transmitted light 

The illuminating light diffuses into the 
substrate and shines through the substrate. The 
system must be able to distinguish the trans- 
mitted light from the defective leakage light. 
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Fig. 30—Structure of roller mask. 


The first two requirements were solved by using 
a roller mask and the third requirement by using 
the self strobing technique. 


6.2.3 Complete masking 

The mask must be able to block the 
illuminating light below the 1077 level. A roller 
mask was developed for the masking. Figure 30 
shows the configuration of the roller mask. The 
roller is a cylinder with a diameter of 2-3 mm. 
The surface is covered with a black light absorb- 
ing material. It contacts the PWB firmly with 
spring action, and rolls when the PWB moves. 
The roller mask blocks the illuminating light 
completely and is capable of following the 
warpage of PWB even as large as 1.5 mm. 

6.2.4 Self strobing 

Leakage light detection sensors cannot 
distinguish between defective leakage light and 
PWB transmitted light. Therefore the self strob- 
ing method was developed to distinguish these 
two types of light. The principle of the method 
is shown in Fig. 31. CCD sensor 1 detects the 
light transmitted directly through the hole. The 
signal is delayed for At by the delay circuit. The 
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X-Y table Illumination X-Y plotter 


Fig. 32—System configuration. 


delay time is set to the time required for the 
PWB to travel from point A to point B. When 
the PWB reaches point B, CCD sensor 2 de- 
tects leakage light from the defective through 
hole. The delayed hole position signal and the 
leakage light signal are input to an AND gate. 
The output of the AND gate indicates the 
presence of a defect. The method can inspect 
only through hole positions automatically, and 
has the advantage that any PWB can be in- 
spected without prior knowledge of the through 
hole positions. 


6.3 System configuration 

The system configuration is shown in 
Fig. 32. The line sensor is made up of 1024 
elements. The PWB travels 10 mm/s on an X-Y 
table. A 300 by 300 mm area can be inspected 
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Substrate 


Normal through hole 


Pad without through hole 


Normal through hole Defective through hole 


Fig. 33—Output of defective through hole. 


Fig. 34—Inspection system. 


in 3 min. The results are shown in Fig. 33. The 
green areas show leakage light signal and the 
white areas show through hole to location signal. 
The two white spots at the left show normal 
through holes, and the black square at the top 
right shows the pad without a through hole. 
The red spot under the right shows a defective 
through hole where the through hole position 
signal and leakage signal coincide. 

The outside view of the inspection system is 
shown in Fig. 34, and the main specifications of 
the system are listed in Table 5. 

The system has been installed at a Fujitsu 
plant. 
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Table 5. Main specification of the through hole 
inspection system 


Items © Specification 


PWB 


7. Mounted components inspection 

Boards with the components mounted on 
PWB are conventionally inspected by using 
electrical function testers and human inspectors. 
The electrical testers can detect defects such as 
incorrect orientation and missing ICs. However, 
a tester cannot detect missing components, 
incorrect orientation, and inadequate insertion 
of capacitors and resistors. The defects do not 
affect functionality, but they decrease the 
reliability of the PCBs. Incorrect orientation of 
capacitors can significantly effect the system, 
that is, the capacitor may break down after 
some usage. The break causes a short circuit and 
can cause the whole system to malfunction. The 
orientation of capacitors must thus be inspected 
accurately. This is accomplished by visual 
inspection, and has been indispensable to assure 
the reliability of PCBs. Visual inspection by 
humans requires much time to achieve the 
necessary reliability. The method cannot adapt 
to the rapid expansion of production quantity. 
An automated inspection method has been 
developed to inspect capacitors. 


7.1 Requirements 

There are two requirements for implement- 
ing this system. 

The first requirement is to test the com- 
ponents with high speed. A new component 
must be inserted by automatic inserting machine 
every 0.5 s. Inspection must be done within this 
time frame. 

The second requirement is high speed image 
processing with a microprocessor to make the 
system cost efficient. Data of a 3D-image is 
much larger than for a 2D-image. For example, 
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Fig. 35—Example of defects. 


a 3D-image with a 256 by 256 pixels resolution, 
and 100 pixels in the vertical direction requires 
4 Mbits, while a 2D-image requires only 64 kbits. 
It is assumed, however, that a microprocessor 
can handle image data less than 0.5 kbits within 
the required 0.5 s. 

High speed posture detection algorithm 
using a microprocessor is required. 


7.2 Sensing method 

7.2.1 Requirements for sensing method 

There are three requirements. 

1) Three dimensional posture measurement of 
black components. 

2) Verification mark printed on the com- 
ponents. 

3) Inspection time of less than 0.5 s per 
component. 

Figure 35 shows the types of several defects 
that must be detected. 

Incorrect orientation and _ position are 
detected by making use of the marks painted on 
components that denote polarity and spec- 
ifications. The marks are usually painted with 
a bright color such as white, and can easily be 
detected with the gray scale images sensed with 
a TV camera. 

Insufficient insertion occurs when the lead 
pin is inserted in the through hole as shown in 
Fig. 35-c). However, gray level images taken by 
TV camera cannot detect the defect. This type 
of defect can only be detected by measuring the 
height of the components. 

Incorrect position can be detected when the 
mark position is seen exceeding the mid position 


between adjacent insertion holes. However, 
components are sometimes inserted at an angle 
tilting as great as 20 degrees. In this case, the 
mark exceeds the limit position. Therefore, 
detecting the only mark position does not 
determine whether the components are inserted 
correctly. When components are inserted in 
adjacent positions and overlap the correct 
position, they may seem correct when viewed 
from above. This means that the components 
posture must be measured to detect wrong 
position. 

Two sensing methods, reading marks painted 
and detecting 3D posture of the components are 
required. 

Conventional 3D posture measurement use 
such techniques as stereo vision and grid pattern 
projection. 

Stereo vision is similar to human 3D 
perception. Two TV cameras are set at different 
positions and are aimed at the same object. 
There are difficulties in detecting black capac- 
itors because little light is reflected. PCBs use 
many black components, so this method is not 
very effective. The processing is also quite slow. 

The pattern projection method uses many 
sheet-type lights and projects them on the 
object. The image obtained shows the structure 
of the object. It takes a long time to reconstruct 
the posture of the object, because there are 
several lines in the image. 

7.2.2 Principle of 3D posture detection 

A new sensing method called the “Two step 
structured light method” has been developed. 
This method can measure 3D posture using the 
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Fig. 36—Principle of Two-step laser light structured 


method. 
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least amount of information without being 
affected by colors. Figure 36 shows the 
configuration of the sensing system. Two laser 
lines are projected on the components with 
a sequentially tilting angle. The laser lines are 
detected with a TV camera. The intensity of the 
laser light is high enough so that even black 
components can be detected. This structured 
light method reduces the number of pixels from 
4 Mbits to 2 kbits. 

The advantage of this method is acquiring 
a 3D posture of the component with the least 
amount of information being required. 

The principle of the method is described as 
follows. 

First, a laser beam strikes position A on the 
components, and the edge position of line /, 
can be detected. After calculating the displace- 
ment and tilting angle of line /,, the 3D 
coordinate of edge points ‘‘a” and ‘‘b” can be 
calculated. 

Next, a second beam is projected to position 
B which is located a known distance from 
position A. The second line determines the 3D 
coordinate of positions ‘“‘c’” and ‘“‘d’’, After 
detecting the 3D coordinate of four positions, 
“a”, ““b”, “‘c’, and “d” of the top surface, the 


orientation and position of the plane can be 


Fig. 37—Laser line image. 
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Fig. 38—Observed image with scattered light. 
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Fig. 39—Configuration of optical system. 


calculated. Figure 37 shows the detected struc- 
tured light on the surface points between a and 
b, and corresponded to the structured line image 
of point a and b shown with a’ and b’. 

The structured light projection method 
enables detection of the 3D posture of com- 
ponents without being affected by color. 
However, it cannot detect marks painted on the 
components. 

7.2.3 Scattered light illumination 

Scattered light illumination is used to detect 
these marks. The components are sometimes 
tilted at angles up to 20 degrees, and the flat 
illumination is needed to illuminate the surface 
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without causing shade. Scattered light is used to 
illuminate the surface. The color of the scattered 
light is blue because it can easily be discrimi- 
nated from the red He-Ne laser light. The 
detected surface is shown in Fig. 38. 

Figure 39 shows the configuration of the 
detection system. The optics for two laser beams 
are arranged to illuminate the components so 
that the lines are perpendicular to each other. 
One of the beams is selected depending on the 
orientation of the components. Two shutters are 
used to select these beams. Each line is detected 
with corresponding TV camera A or B. The 
mark image is sensed by a TV camera C 
mounted above the components. 


7.3 Inspection algorithms 
The algorithms detect component orien- 
tation, mark position, and determine whether 
the components are correctly inserted. 
There is an algorithm for posture measure- 
ment and one for orientation measurement. 
7.3.1 Posture measurement 
Posture can be calculated from the two step 
structured light method. This method has the 
following advantages. 
1) Processing time of less than 0.5 s/component 
2) Segmentation of the line image on the 
components from the surrounding com- 
ponents images 
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Segmented area 


Fig. 40—Segmented line. 


3) Wide dynamic range 

Reflected light intensity can vary as much as 
ten times depending on the optical character- 
istics of the surface such as color and roughness. 
The inspection algorithms should not be 
affected by this. 

To meet the first requirement, a window is 
set at the line image of the laser light to limit the 
detection area. The window position can be 
calculated from the size of the components 
being inspected that are known prior to the 
inspection. The window reduces the image data 
from four million pixels down to two thousand 
pixels to shorten the processing time to 0.5 s. 

For the second requirement, the continuity 
of the line is calculated to segment the line 
image. The continuity is detected by calculating 
the angle of tilt and the length of the same angle 
region. The length of each segment is compared 
with the tolerance. Only the one segment 
nearest the correct one is selected as the image 
of the component that is inspected. Figure 40 
shows the segmented results of the line image. 

The third requirement is satisfied by 
detecting the brightest points of the line pattern. 
The line intensity is checked perpendicular to its 
direction and the center of the line is selected. 
The peak detecting method widens the dynamic 
range to ten times of that without it. This 
process works well even when the line intensity 
exceeds the limit of the sensor’s dynamic range. 
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Fig. 41—Double window. 
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Fig. 42—System configuration. 


7.3.2 Polarity detection 

This algorithm detects whether the marks 
are seen at the correct positions. This method 
checks the mark orientations using two variable 
windows. One window is set where the correct 
position mark may appear and the other window 
is set where the wrong mark will appear. The 
window positions are displaced to the angle of 
tilt as measured by the structured light method. 
In each window, areas brighter than the 
threshold are detected and compared between 
the two windows. When the bright area in the 
correct window is wider than the wrong one, 
then the orientation is judged to be correct 
because the mark area is usually wider than that 
of noise patterns. Figure 41 shows an example 
of a detected image in the windows. 


7.4 System configuration 
Figure 42 shows the system configuration. 
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Fig. 43 Experimental results. 
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Fig. 44—Inspection system. 


The images that are taken from the three TV 
cameras are stored in frame memories. The 
frame size is 256 by 256 pixels with 16 gray 
levels. A frame is taken in 33 ms. An 8088 
microprocessor is used for calculating the 
component orientation, and the angle of tilt. 
An Intel 8087 numeric coprocessor is used for 
triangular calculation. 


7.5 Inspection characteristics 

Figure 43 shows the accuracy of calculating 
the angle of tilt. These results were obtained 
after 20 measurements. They show the mean 
and standard deviation, and the dotted lines are 
at the 30 level. Samples are tilted from —25 
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Table 6. Main specification of the component 
inspection system 


Item Specifications 


Components Rectangular capacitors 


Size of PCB 300 x 400 mm, ¢ = 1.2-2.5 mm 
: Horizontal: 0.1 mm, 
ie atin Vertical: 0.2mm 
Missing, incorrect position, 
Defects inadequate insertion, 


incorrect orientation 


Inspection speed 0.5 s/component 


Size LSX1S 612 xb xh) 


degrees to 25 degrees. The required accuracy is 
+10 degrees. Figure 43 shows that accuracy 
requirement is satisfied. There are some errors 
because the shape of the components are not 
rectangular cubics. 

A practical inspection system has been 
constructed as shown in Fig. 44. The missing 
rate is zero percent, and the detection rate is 
99.7 percent. The inspection time is 0.5 s per 
component. The system is now on the produc- 
tion line at a Fujitsu plant. 

The main specifications of the system are 
listed in Table 6. 


8. Conclusion 

The technologies for automatic PWB and 
PCB visual inspection have been developed. 
Based on these technologies, automatic inspec- 
tion systems had been developed by Fujitsu, and 
working systems have been installed at the 
Nagano, Akashi, and Oyama plants since 1979. 
These systems have been effective in achieving 
high reliability of PWBs and PCBs which are 
used in large scale computers, etc. 

Up to now, further developments in auto- 
matic inspection technologies have been con- 
tinued until the inspection speed increased to 
ten times, and the resolution as fine as 5 ym. 
The technologies are being implemented in the 
latest inner layer pattern inspection system. 

In the future, PWBs will have finer patterns 
on a larger scale, and will demand further high 
reliability. Inspection technologies will play 
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This paper describes the technical problems encountered in making a solid CAD system fit 


for practical use, and the solutions to these problems. The paper is based on the results of 
the development of a three-dimensional solid CAD system, |CAD/SOLID. 
The following three techniques and functions were developed to make a solid CAD system 


practical: A man-machine interface that allows a designer to create and operate solid models 


using color shading images, functions that allow two-dimensional drawing data to be used 
directly for three-dimensional shape operations, and functions to develop and add specific 


functions using a solid model. 


These functions do not require a knowledge of the internal data structure of a solid model 


or geometrical calculations. 


1. Introduction 

Solid modeling techniques first appeared 
more than ten years ago. These techniques 
have been widely known to be both effective 
and useful for general-purpose applications. 
But solid CAD systems, which are CAD systems 
that use the solid modeling technique, have not 
become popular?*?). Solid CAD systems are 
presently regarded as systems to be used only 
for a specific business by a specific user. The 
aim of this paper is to illustrate the technical 
problems that have prevented solid CAD systems 
from becoming popular, and to discuss the 
solutions to these problems. 

Research on solid CAD systems usually 
means research on solid modeling techniques. 
Examples are the interference detection between 
solid models, and control of the shape of a 
sculptured surface. These are very important 
areas of research, and are _ necessary for 
developing a powerful solid CAD = system. 
But a powerful solid modeling technique does 
not always result in a practical solid CAD 
system. There are other important areas of 
research that tend to be overlooked, for 
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example, research on the man-machine interface 
for interactive shape operations. It is important 
for a designer to create and operate on solid 
models easily and quickly. A powerful solid 
modeling technique and an easily used man- 
machine interface allow a_ practicable solid 
CAD system to be developed. This paper 
discusses the second problem, the design of a 
natural and easy-to-use man-machine interface 
for a practical solid CAD system. 

Our research started with a trial and 
evaluation of solid CAD systems, which clearly 
exposed the problems, and are discussed in 
chapter 2. Three new techniques were developed 
from the analysis of the problems, which are 
described in chapter3. The techniques were 
applied to the three-dimensional solid CAD 
system, ICAD/SOLID, whose features are 
described in chapter4. The application of 
these techniques to the design of a telephone is 
described in chapter 5. 

The conclusions are given in chapter 6. 
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Fig. 1—A wireframe picture. 


2. Technical problems with the solid CAD sys- 
tem 
Tests and evaluations of solid CAD systems 
have exposed the three main problems that 
make solid CAD systems impractical. 


2.1 A method of interaction 

To design products effectively, it is 
important for a designer to be able to create 
and operate solid models easily and quickly. 
Easy and quick operation of an interactive solid 
CAD system depends on the form of the display 
and the method of interaction. There are three 
main considerations for the man-machine 
interface: the display method, the element 
indication method and the delay in displaying 
the image. 

2.1.1 Display method 

Usually, a man-machine interface based on 
wireframe pictures is used for an interactive 
shape operation, because such images require 
fewer calculations. The advantage of wireframe 
pictures is therefore little delay in displaying 
the image. An example of a wireframe picture 
is shown in Fig. 1. However, as a designer 
geneates and operates solids in a solid CAD 
system, and shape operation in a solid CAD 
system is based on edges, faces and solids, 
it is difficult for the designer to understand 
solid shapes from wireframe pictures and 
Operate upon them. 

Color shading makes it easy for a designer 
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a) Crankshaft 


b) Printed circuit board 


Fig. 2—Color shading images. 


to understand solid shapes. Examples of color 
shading images are shown in Fig. 2. A man- 
machine interface based on images with color 
shading is suitable for solid shape operation. 
But as many calculations are needed to generate 
color shading images from solid models, the 
man-machine interfaces based on color shading 
images are seldom used for performing solid 
shape operations. The technique is normally 
used only for verifying finished shapes. 

The above method, an interactive shape 
operation in wireframe pictures and a verifi- 
cation of the finished shapes using color shading 
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Fig. 3—A face indication by hits of two edges. 


images, is the usual one. However, it makes the 
man-machine interface of a solid CAD system 
difficult to use. 

2.1.2 Element indication method 

A designer must operate upon and indicate 
edges, faces, and solids in a solid CAD system. 
The element indication method depends on 
the display method. In wireframe pictures, 
a designer can indicate an edge directly in one 
hit action. However, a face cannot be indicated 
directly, as the faces are not actually displayed. 
It therefore takes two actions to indicate a face, 
for example, two of the edges on the face must 
be hit (see Fig. 3). In color-shaded images, a face 
can be indicated in one action, while an edge 
needs two actions. 

For simple and quick operation, both an 
edge and a face should be indicated in one 
action. But neither wireframe pictures nor color- 
shaded images can satisfy this requirement. 

2.1.3 Delay in displaying the image 

Little delay in displaying the image is 
important in interactive solid CAD systems. 
When a shape operation has been performed, 
a designer must confirm the result of the 
operation in real time. The designer evaluates 
the shape operation by looking at the updated 
image, and decides upon subsequent action. 
In order not to interrupt this stream of thought 
and action, the updated images should be 
displayed in real time. 
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Fig. 4—A boolean operation. 


As described in subsection 2.1.1, wireframe 
pictures can nearly be displayed in real time. 
But it takes a lot of practice for the designer to 
understand solid shapes from wireframe 
pictures. And although it does not take long to 
understand solid shapes from  color-shaded 
images, such images require many calculations 
and cannot be displayed in real time. 


2.2 Shape operation method 

The natural and efficient solid CAD system 
depends on the method used for the shape 
operation. The ideal operating method matches 
a designer’s natural way of thinking while 
an unnatural shape operation method will 
make an interactive solid CAD system im- 
practical. 

The starting point for designs is still the 
drawing. Many designers still prepare their plans 
on drawing boards. They do not design three- 
dimensional products straight away. They put 
their ideas into drawings, and then’ they 
reconstruct solid images in their minds from 
these drawings. 

The shape operation functions of a solid 
CAD system usually consist of a direct gener- 
ation of solid primitives and a boolean operation 
between two primitives. An example of a 
boolean operation is shown in Fig. 4. Three- 
dimensional products are treated directly as 
solid primitives in a solid CAD system. But, 
this standard method does not always fit the 
designer’s natural way of thinking. This is 
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Fig. 5—Examples of input faces. 


the second problem that often makes solid CAD 
systems impractical. 


2.3 The method of adding specific functions 

For a solid CAD system to be effective, 
the user must be able to easily add user-defined 
functions to the system. It is desirable that users 
develop these specific functions and add them 
to their systems themselves, as they understand 
their own business best. 

However, as a lot of technical knowledge 
is needed to develop the programs, it is difficult 
for a user to develop specific function programs 
using a solid model. Knowledge of winged-edge 
data structures and boundary representations 
is needed, for example. This is the third problem 
that makes present solid CAD systems im- 
practical. 


3. Techniques for solving the problems 
The following techniques were developed 
for solving the problems described in chapter 2. 


3.1 Techniques for interaction 

As described in section 2.1, man-machine 
interfaces based on color-shaded images are 
desirable for solid shape operation. However 
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for easy and quick operation, edges and faces 
should be displayed in the same picture. It is 
highly desirable that the solid-model display 
used for the shape operation overlaps the wire- 
frame pictures that represent edges, with color- 
shaded images representing the faces. A three 
dimensional graphic display with a function 
for eliminating hidden surfaces is necessary 
for displaying overlapping pictures in real 
time. The functions described below have 
been developed. It is assumed that a three- 
dimensional graphic display with a hidden 
surface elimination function is used. Solid 
model data are converted into polygon data 
that the display can process. The polygon 
data are transmitted from the host computer 
to the display, which eliminates hidden surfaces, 
and displays the solid model. 


3.1.1 Polygon division method 

A polygon division method can transform 
a solid model given in boundary representations 
to three-dimensional polygon data form. It 
divides curved faces that include holes into 
planar polygons without holes, and curved 
edges on the boundary of the face into a chain 
of straight lines. 
1) Algorithm 

i) Input 

Input data consists of the face, edges and 

vertices of the boundary representations. 

This method can divide a planar surface, 

a quadratic surface, or a sculptured surface, 

and faces can have holes. Examples of 

input faces are shown in Fig. 5. 


ii) Division of edges 

Curved edges on the boundary of a face 
can be divided into straight lines. For 
example, a circle edge can be divided into 
a number of short lines, as shown in Fig. 6. 
The width of each division can be controlled 
according to which the designer thinks is 
more important: a rapid display, or a precise 
display. This chain of straight lines is used 
for wireframe pictures and it becomes the 
input of the step “division of a face”’. 
iii) Transformation into (u, v) space 

Short straight lines in three-dimensional 
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Fig. 6—Division of an edge. 


real space (x, vy, z) can be converted into 
two-dimensional parametric space (u, V). 
This step transforms a curved face in (x, y, 
z) space into a polygon in (u, v) space, as 
shown in Fig. 7 and enables all the faces to 
be processed in the same way. In the next 
step, the polygon is divided into a number of 
triangles and rectangles that the display can 
process. 
iv) Division of a face 

It is important to fit wireframe pictures 
to the boundaries of the faces in overlapping 
pictures. If wireframe pictures do not fit 
exactly onto the boundary of the face, they 
may be covered by the face, and be put out 
of sight. Likewise, if the boundary of one 
face does not exactly fit the boundary of 
a neighboring face, there will be gaps be- 
tween the two faces. Division of a face is 
usually done by a mesh method, which 
divides faces into grid squares by respective 
meshes of the faces, as shown in Fig. 8, 
repeats division until polygons in all meshes 
become triangles or rectangles. This means 
that the boundary of a face does not always 
fit the boundary of a neighboring face. To 
solve these problems, this algorithm gener- 
ates color shading data for the neighboring 
faces from only the wireframe data of the 
common edge. Triangles and rectangles are 
generated in this step. 
a) Generation of a triangle from an acute 

vertex 


Fig. 7—Transformation from (x,y,z) to(u,v). 


Fig. 8—Division by the mesh method. 
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Fig. 9—Generation of triangles. 


Check the interior angle of every vertex 
of the input polygon. If the interior angle of 
the vertex is acute, generate a triangle 
polygon from the vertex and neighboring 
two vertices, as shown in Fig. 9. Delete the 
vertex from the chain and if there is no 
acute-angled vertex, go to the next step. 

b) Generation of a right-angled triangle 
from a hypotenuse 

Check the angle of inclination of every 
edge of the input polygon. Generate a right- 
angled triangle polygon from the edge. 
If the edge is not parallel to the u-axis and 
y-axis, a new vertiex is generated in the input 
polygon, as shown in Fig. 10. Add the vertex 
to the chain. The edge is the hypotenuse of 
the right-angled triangle, and the generated 
vertex is at the right-angled corner. If all 
edges are parallel to the u-axis or the y-axis, 
go to the next step. 

c) Generation of a rectangle 

Arrange the vertices in increasing order 
of yv-coordinates. Arrange the vertices that 
have same y-coordinates in increasing order 
of u-coordinates. 

Search this chain of vertices, cut at the 
same  -coordinates, generate rectangle 
polygons and delete the vertices from the 
chain, as shown in Fig. 11, if there is no 
vertex in the chain, go to the next step. 
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Fig. 10—Generation of right-angled triangles. 
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Fig. 11—Generation of rectangles. 


d) Subdivision of a rectangular polygon 
Subdivide a rectangular polygon at even 

intervals of u-coordinates into smaller 

rectangular polygons, as shown in Fig. 12. 

v) Transformation into (x, y, Zz) space 
Transform polygon in (u, v) space into 
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Fig. 12—Subdivision of rectangle polygons. 


(x,y, 2) space. 

Add normal vectors to the vertices of the 
polygons. 

vi) Output 

The output data consist of polygons and 
polygonal lines. Polygon data consist of 
triangles and rectangles that have no holes. 

The vertex of the polygon has a normal 

vector, and the vector is used for color 

shading. Polygonal line data exactly fit the 
boundary of the faces that are approximated 
to planar polygons. 

3.1.2 Extraction of updated faces 

It is important in order to have rapid display 
to minimize the amount of display data that 
must be transmitted from the host computer to 
the graphic terminal. It is not necessary to 
transmit display data for a whole solid model 
when a shape operation is performed. Only the 
display data for the faces that are transformed 
A or generated B- should be transmitted, 
as shown in Fig. 13. 

A function for extracting the updated faces 
has been developed, which is built into the basic 
subroutines constituting the shape operation. 
During a shape operation, these subroutines 
write the records of generated, deleted and 
transformed elements into a table. The faces 
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Fig. 13—Extraction of transformed A_ or generated 
B_ faces. 


that must be transmitted are selected according 
to this table. 


3.2 Technique for natural operation 

According to the example of a boolean 
operation between simple primitives shown 
in Fig.4, it seems that a boolean operation 
is easy to use. But some designers who have 
tried the solid CAD systems have pointed 
out that it is very difficult to generate cutter 
solid models for a boolean operation. A cutter 
solid model is necessary in a boolean operation 
to indicate the shape to be deleted. However 
it is unnatural and difficult for a designer to 
generate a complex cutter solid model directly 
without using drawings, although it is effective 
to generate a cutter model using drawing data. 
And it is more effective to use drawing data 
directly for a shape operation. 

Addition of the data of the third coordinate 
is necessary in order to use two-dimensional 
drawing data for a three-dimensional shape 
operation. 
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Conversation processing 


Specific functions 


Standard functions 


Display processing 


Terminal function 


Hidden-surface processing 


Solid model 


Rotation, scaling, translation 
Shading 
4096 colors 


DBMS 


MAC: Model ACcess Routine SPL: Shape Processing Library GCL: Geometrical Calculation Library 


Fig. 14—System configuration of ICAD/SOLID. 


Although it is difficult to add the data for 
the third dimension automatically, it is possible 
to add this data using the depth information of 
the viewport in which the drawing is displayed. 
Moreover, the drawing data is stored in the 
same format record as the three-dimensional 
edge for a quick transformation into a solid 
model. 


3.3 Techniques for adding specific functions 

It is difficult for a user to develop specific 
functions using a solid model. Therefore, the 
subroutines that are open to a user should 
not require detailed knowledge of the internal 
data structure of a solid model and geometrical 
calculations. 

The subroutines were developed, based 
only on data structures that a user can under- 
stand through overlapping pictures. The data 
structure consists of solids, faces, edges, and 
vertices. The subroutines combine basic sub- 
routines that use the internal data structure, 
and put the internal data structure out of sight. 


4. Features of ICAD/SOLID 


The three-dimensional solid CAD system 
ICAD/SOLID is a new part of Fujitsu’s 
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integrated computer aided design and manu- 
facturing system (ICAD). Using the techniques 
described in chapter 3, the following functions 
have been realized in ICAD/SOLID. The system 
configuration of ICAD/SOLID is shown in 
Fig. 14. 


4.1 An advanced man-machine interface 

All the operations of ICAD/SOLID are 
carried out by combining color-shaded images 
and wireframe pictures. A solid, a face, and 
an edge can be indicated in a single action, 
consisting of one hit with the pointing device. 
If necessary, the display can be converted 
into color-shaded images or wireframe pictures, 
as shown in Fig. 15. 

The updated image after a shape operation 
is displayed in a few seconds. Using the func- 
tions of a three-dimensional graphic terminal, 
cross sections and internal views of solid models, 
as shown in Fig. 16, can be obtained within 
a second. 

ICAD/SOLID is an advanced man-machine 
interface for three-dimensional shape operations 
which is compatible with a designer’s thought 
and work. 
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Fig. 15—A change of a display type. 


a) Internal view of telephone 


b) Cross section of telephone case 


Fig. 16—Cross section and internal view. 
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Cutter lines 


Fig. 17—MILL function. 


4.2 Shape operation based on drawings 

ICAD/SOLID can receive drawing data 
from ICAD/SDS” (integrated computer aided 
design and manufacturing system/standard 
design system). Drawing data is used for shape 
operations and generating solid models while 
the MILL function, shown in Fig. 17, is the 
typical function for a shape operation. This 
function adds the depth or height information 
to the drawing data and uses the drawing data 
as cutter lines which indicate the shape to be 
deleted. It does not use a cutter solid model, 
but uses drawings directly. 

The method fits a designer’s natural thought. 


4.3 Structure open for adding user-defined func- 
tions 
A user can incorporate programs written 
in FORTRAN77 into ICAD/SOLID. The user’s 
additional programs can call up standard 
programs from ICAD/SOLID. ICAD/SOLID 
has the following standard subroutine libraries. 
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Outer case 


Allowable mounting height 
of electronic parts 


PC-board 


Fig. 20—Concept of packaging space calculation function. 


Fig. 18—A solid model of Fujitsu telephone GK. 


Investigation into interference and position of components to 
be mounted (hook switch, volume control, etc.) 


u Fig. 21—Output of packaging space calculation function. 


Projected drawings including cross-sections for detailed 5. Application in layout design 
i ncaa aaa ICAD/SOLID has been used for developing 


the layout and structural design of a telephone”. 
A drawing system Solid models of the FUJITSU telephone GK are 
shown in Fig. 18. The standard workflow is 

Fig. 19—Standard workflow. shown in Fig. 19. 


The user department developed the ad- 
ditional programs for the specific functions 
2) Geometrical calculation library with the help of the libraries belonging to 
3) Model access library ICAD/SOLID, which is the parckaging space 
4) Display and interaction library calculation function for printed circuit boards. 

The function divides the surface of a printed 
circuit board into square grids, and then 


1) Shape processing library 
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calculates and outputs the allowable mounting 
height of the electronic parts that are mounted 
on the printed circuit board. The concept 
of this function is shown in Fig. 20, and an 
example of the output is shown in Fig. 21. 

Packaging space is normally calculated by 
hand and the information produced is ap- 
proximate. However using ICAD/SOLID and 
this additional function, the packaging space 
can be calculated accurately. 


6. Conclusion 

With the development of a man-machine 
interface for three-dimensional shape operation, 
a designer can design a product of complex 
shape with high design quality more efficiently. 

The new shape operation functions have 
made it possible for a designer to use drawing 
data directly for solid modeling. 

A program interface for adding specific 
functions has made it possible for the user 
department to easily develop application 
programs using a solid model. Having discussed 
the man-machine interface for solid shape 
operations, the next step is to develop the man- 
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machine interface that raises the efficiency 
of the application of solid models. 
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One of the inherent disadvantages of the Hough transform, its long computation time, is 
preventing this transform from being used in real-time applications. This paper presents 
an architecture that was developed to address this problem and implemented in an experi- 
mental hardware model for real-time straight-line detection. 

This paper also provides examples of applying this architecture to the automatic inspection 
and measurement of objects in factories and laboratories. In these examples, it took less than 
one second to extract straight-line parameters from captured pictures. 


The highly parallel architecture proposed here for the Hough transform does not require 


any multiplication or trigonometric calculations at run time. 


1. Introduction 

The Hough transform for recognition can ef- 
fectively process noisy backgrounds and bound- 
ary gaps. It was first applied to the detection 
of straight lines!) and has since been generalized 
to successfully detect parametric curves such as 
circles, parabolas, and ellipses. Recent literature 
about the Hough transform includes many 
works related to applications such as recognizing 
nonparametric curves and 3-dimensional objects, 
determining the orientation of rigid bodies, and 
obtaining motion parameters for moving ob- 
jects?) 3), 

In spite of the excellent recognition ability 
of the Hough transform, such disadvantages as 
long computation times and large memory 
requirements prevent its being used in real-time 
applications??. In general Hough transforms, 
computation time and memory requirements 
increase exponentially with the number of 
dimensions in the parameter space*” ©). 

Silberberg reports the execution time of 


a Hough transform program for straight-line 
detection on a VAX 11/7857). Only integer 
arithmetic was used with cosine and sine values 
stored in a table. The program he ran took 
about seven seconds to process a 521 x 512 
image without peak detection. 

One of the objectives of this paper is to 
present an architecture developed to address the 
time and memory load problems. This architec- 
ture has been implemented on an experimental 
hardware model for real-time straight-line detec- 
tion. 

The underlying premise of our study is that 
large memory requirements are not serious 
handicaps at current levels of memory integra- 
tion technology for straight-line detection in 
two-dimensional arrays, which is the most 
fundamental and widely applied use of the 
Hough transform. 

As can be seen from the literature, major 
applications of the Hough transform are still in 
the field of straight-line detection. Some of 


This paper has been revised from the article “High-speed Hough transform processor and its applications to 
automatic inspection and measurement” as it appeared in the proceedings of the IEEE International Conference on 


Robotics and Automation, April 7-10, 1986. 
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these applications include gauge inspection®?, 


recognition of handprinted Hebrew characters”), 
inspection of VLSI resist patterns!®), detection 


11) and vehicle guidance 
12)-15) 


of seismic patterns 
through road edge detection 

The time required for our model to make 
a histogram in a 512 x 512 (p, 6) space for 1 024 
feature points is 0.26 s. These results can be 
compared to those of the report”) after correc- 
tion by appropriate extrapolation. Execution 
time is linear for the number of the feature 
points and partitions in the @ coordinates. The 
report assumes 3 percent feature points in the 
512x512 image and five degrees for @ reso- 
lution, which implies 12.9 s for our conditions. 

Chapter 2 describes the processing flow of 
the Hough transform, and chapter 3 presents our 
design objectives for the hardware implementa- 
tion of Hough transforms. Chapters 4 and 5 
describe the optimum parameter space reso- 
lution and quantization errors. Chapter 6 shows 
the system configuration of our experimental 
hardware model, chapter 7 shows that the 
complementary Hough transforms are suitable 
for the developed processor architecture, and 
chapter 8 proposes a parallel architecture for 
Hough transforms. Chapter 9 presents examples 
of applying this architecture to automatic 
inspection and measurement tasks. 


2. Processing flow 

This chapter first defines the symbols used 
to set the coordinate systems of the image and 
parameter space. It then describes the param- 
eterization of straight lines. 

Raster-scanned TV signals are sampled to 
assemble an image space consisting of Nx N 
pixels. Assuming an even aspect ratio, this space 
is a square. 

A pair of horizontal and vertical addresses 
(H, V) is often used to indicate a location in a 
raster-scanned picture. This pair of addresses is 
referred to as the raster address. We use the 
following normal parameterization proposed by 
Duda and Hart!” to specify the straight line: 


p=x;cos@+ y; siné, 
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where, x; and y; are the relative values from the 
center of the picture at the i-th feature point in 
the picture. 

Assuming the raster address of this feature 
point to be (H;, V;), the address values are given 
by the relations 


N 
ss ees es ee 


where N is assumed to be an even number. In 
our system, N = 256. p is the distance of the line 


ae : N N 
from the origin with raster address Vara) 


and @ is the angle between the normal and the 
xX axes. 

The sinusoidal curve (1) in the parameter 
space represents all straight lines passing through 
the i-th feature point. 

After all feature point curves (1) have been 
drawn in the parameter space (p, 6), the curves 
corresponding to collinear points in the picture 
should intersect at one point. Assuming the 
coordinate of this point to be (pp, 4), the line 
that passes through the collinear points can be 
represented in the image space as follows: 

Pp =X COBO, +VSINGD,, = vanes (2) 

This reduces the problem of detecting the 
straight lines in the image space to a problem of 
searching for the intersections of sinusoidal 
curves in the parameter space. 

The method used to search for these inter- 
sections is to quantize the (p,6@) space and treat 
it as a two-dimensional accumulator array. First, 
the array is cleared. Then, all cells of the array 
along the curves (1) are incremented by one 
unit. Once sinusoidal curves are accumulated 
for all feature points, the intersections appear 
as peaks. 


3. Design approach 
3.1 High resolution 

Coarser quantization shortens processing 
time and reduces the demand for accumulator 
array memory space. Furthermore, the resulting 
histogram of the locus of the sinusoidal curve 
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will show legible peaks even for boundaries or 
edges of objects with poor linearity. The dis- 
advantage of coarser quantization is that there 
are significant errors in parameters obtained 
from boundaries or edges with good linearity. 
We have adopted a design for measuring 
objects with good linearity using high resolution 
by selecting the optimal step size. 
Two-dimensional arrays are sufficient for 
detecting straight lines, and memory utilization 
requirements are not a serious problem with 
current memory integration technologies. The 
same histogram that results from a transform 
with a coarse step size can be obtained by merg- 
ing multiple cells in postprocessing after 


histogram creation. This merge processing is 
a simple local operation. Our concern was 


that it isimpossible to obtain fine histogram 
structures from histograms with coarse step 
sizes. 

Over-sampling of the (p, 6) space, however, 
must be avoided. Research is now being done to 
find the optimum relationship between A@ and 
Ap'”. Chapter 4 introduces similar work 
resulting from a different point of view. 


3.2 Parallel operations 

We have separated sine curve generation and 
histogram generation so that they can be done 
in parallel. 

Calculating the increment for a single feature 
point cannot begin until p, has been calculated 
for that point. Normally, p calculation and 
increment processing cannot be done in parallel. 
They can be executed in parallel, however, if the 
feature points are input consecutively. p, for 
the i+ 1-th point can then be calculated while 
the i-th point increment processing is being done 
because these two operations do not share hard- 
ware resources. 


3.3 High-speed data transfer 

In the edge extraction section, convolution, 
thinning, local maxima detection in nearby 
pixels and zero-crossing detection can be per- 
formed at video rates!®). Because the Hough 
transform section cannot complete processing 
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for each address within one sampling period, the 
addresses are temporarily buffered by a FIFO 
(first-in first-out) buffer. The histograms are 
transferred back to the space differential pro- 
cessing unit for edge extraction and _ local 
maxima are searched for at a video rate. 

The sets of (p, 8) coordinates detected by 
hardware are also sent to the FIFO buffer in 
the order they are found along the raster scan. 
This means that the FIFO buffer is used twice 
during the complete process. The type of data 
passed through the buffer changes. It first carries 
addresses in the image space and then carries 
addresses in the parameter space. 


4. Parameter space resolution 

This chapter describes the optimum relation- 
ship between A@ and Ap. 

If the @ pitch is large in the discrete param- 
eter space, p values must sometimes be inter- 
polated to draw a sine curve. Interpolation is 
required when the difference in p values at 
adjacent @ coordinates exceeds Ap. This is 
likely to occur when the curve crosses the 0 
axis and x; and y; are both large. Hardware 
interpolation has the disadvantage of requiring 
either additional circuits or more complex 
processing, both of which lower processing 
speeds. It is possible, however, to find a Ap 
value that makes interpolation unnecessary. 

To achieve this, even when 


dif (p)=p (On+1 5 Mis Va) — p (Ons Mie Vi)s 


is maximized requires that the following equa- 
tion be satisfied. 
0p TT / xe +y;? 


dif _ = e 
(p) An Mo 


Tn 
cos (—— +6 9)< Ap, 
Ms 0) 
for all x;, yj; and m, where we approximate the 
difference in ¢ by differentiation and n is con- 
sidered continuous. 


s 2 a. 
Because the maximum value of \/ x;7+ y;’ is 
NAx 


Vt 


, Equation is simplified as follows: 
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TmNAXx < Ap (4) 

No JD Tg SSO, Me saieecrans eterecer 

The optimum relationship between N, and 
Ap is where equality holds. 

So far we have assumed that cos 6, and 
sin @, are of infinite precision. Now, we must 
check their accuracy. If the rounding error of 
cos 0, is 6C and that of sin 6, is 6S, the square 
of the p value rounding error satisfies the follow- 
ing Schwartz inequality: 


(5p)? =(x; 6C +9; 6S)? < 
(x;? +y;?) (8C? + 8S? )<.N? (Ax)? 8?, 


where 5? is the maximum value of 6C? or 
5S?. For 5p to be less than Ap along any curves, 


5 must be less than 


x 
Finally, the maximum number for each 
cell of the accumulator array can be determined. 


TV : $16 : D/A converter 
camera A/D con- — 
verter 


16 
eal 


64-Kbyte frame 
memory 


Internal bus 


Edge extraction Edge flag 


section 


Video bus 


The maximum count equals the maximum 
number of feature points that can be contained 
in a bar-shaped window!” whose width is Ap. 
The window contains the maximum number of 
points when it is oriented diagonally in the 
image space. The area of the common portion of 
the image and the window is approximately 
V2 NAx:Ap. There is one feature point in 
each square whose size is Ax, so the average 
Ap 

Ax 
The bit length of the counter that repre- 
sents one cell and must count up to 


A 
jogs afd N= bits. 


number of points in the window is./ 2 .N 


5. Hardware specifications 

This chapter describes how the quantization 
parameters, whose optimum relationships were 
described in the preceding chapter, were reflect- 
ed in the experimenta! hardware model. 


Console 
Monitor 


UNIX 


IEEE 796 bus 


Parallel I/O Control/status 
MC68000 RAM, ROM 


Hough transform section 1 

: j ; | 
Histogram Histogram processing | 
memory section | 
256 Kbytes H(n, pn) = H(n, pn) +1 | 
| 

| 

| 

cos @sin @ table po=%X, cosé@+ ¥ sin 8 

| 

| 

é | 

Buffer 1 "3 p calculation Buffer 2 | 
: section 

| 

eae a ee ee axl 


Fig. 1—System configuration. 
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The size of the image space is 256 x 256, so 
N= 256 and the ratio of Ap to Ax is set to 1. 
In other words, adjacent horizontal or vertical 
lines can be differentiated. Under these con- 
ditions, the number of p partitions can reach 
25644) 2. According to the optimum relation- 
ship between Ap and AO, the @ partition number 

256-7 

i 
requirements, the number of p and @ partitions 
were both set to 512. This results in some 
redundancy because (512 — 256-/2) Nog cells 
will never be used. The step size of @ is coarser 
than the optimum value, so interpolation is 
desirable, but we did not use interpolation 
because Equation (3) is satisfied for most 
of the feature points (i.e., those in a circular 


is To simplify hardware manufacturing 


range within a radius of 163 picture elements). 
Although log) Y2N = 8.5 bits are required 
for each cell of the accumulator array, only 
eight bits were taken. This was done for two 
reasons: first, to ensure compatibility with 
other parts of the image processors or micro- 
processors in the experimental model and 
second, because continuous straight diagonal 
lines are rare. The values cos 6, and sin 6, were 
quantized at eleven bits including signs. 


6. System configuration 

Figure 1 shows the system configuration, 
which comprises three main sections: the edge 
extraction section for extracting edge points 
from TV signals, the Hough transform section 
where sine curve histograms are created from 
extracted edge points according to Equation (1), 
and the microprocessor for general control of 
the total system and its peripheral devices. 


6.1 Edge extraction section 

TV signals are sampled at 6 MHz and quan- 
tized in 8-bit samples. The edge extraction sec- 
tion processes the quantized signals at the video 
rate. This section also has two convolver stages. 
By changing stage combinations, either Sobel 
operation + fine-line binary coding or smooth- 
ing + Laplacian + zero crossing detection can be 
selected (see Fig. 2). The mask is either 3 x 3 
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Local-bus 


Thresholding 


Edge 
flag 


Hake. | 
Be IP 
Hew) > 0 
lJ > 
ep 


Ridge/peak extraction 
Ds 
Fs 
Bs 
=P =p —D; 
Cs Z 
e— > 
As 
—— 


Zero crossing detection 


= 


CM: Comparator 
DS : Data selector 


Fig. 2—Edge Extracting Structure. 


or 4x4 (see Table 1). The scan addresses of 
edge data obtained here are transferred to the 
Hough transform section at the video rate 
through a bipolar FIFO buffer with a capacity 
of 1 kW (1 W= 16 bits). 


6.2 Hough transform section 

The Hough transform section consists of two 
subsections. The p calculation section uses 
Equation (1) to calculate p at each edge point 
and make a 512-word (n, p,) table for each 
edge point. These tables are subsequently used 
by the histogram processing section. The con- 
tents of the histogram memory location ad- 
dressed by (n, p,) are incremented by one for 
n=0 to 511. This procedure is repeated for all 
edge points. The p calculation section consists of 
bipolar ROMs, where sine and cosine values are 
stored to an accuracy of eleven bits, two 50 ns 
multipliers, and a 16-bit adder (see Fig. 3). 


39 


K. Hanahara et al.: A Real-Time Hough Transform Processor 


Table 1. Convolver coefficients 


Fil Coefficients 
ilter 
: 3XK3 4x 4(1) 4x 4(2) 
Smoothing 
Laplacian 
Enhancement 
Gradient 
(Sobel) 
1 
2 4 4 
—2 -4 -4 
—1 -—2 -2 -1 


The histogram memory has a capacity of 
512x512x8 bits. It consists of 64 K-bit 
CMOS RAMs with an access time of 150 ns. 


7. The complementary Hough transform 
Brown proposes a new scheme for Hough 
transforms which uses negative votes as well as 
the usual positive votes”. This complementary 
Hough transform, as he calls it, can eliminate 
off-peak positive votes in the parameter space, 
sharpen peaks, and reduce sidelobe bias and 
variance. For detecting circles, Brown uses 
feature point spread functions consisting of a 
ring of Is surrounded both inside and out by a 
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ring of == 5: 


The complementary Hough transform can be 
considered almost equivalent to performing 
the Laplacian on the final histogram based on 
positive votes. When there is only one feature 
point in the image space, the sinusoidal curves 
consisting of positive and negative votes are 
almost the same as the result obtained by apply- 
ing the Laplacian to the curve. 

Further, because the Laplacian is a linear 
operation and the superposition principle applies 
to the process of making histograms, the 
sequence of these two operations can be re- 
versed freely. 
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aitis Galea Modified complementary Hough transforms 
—— cos 6, : based on the Laplacian are suitable for the proc- 


table 
512 words 


essor architecture. All that is required is to 
perform the Laplacian on the histogram in frame 
memory before performing peak detection. The 
extra processing time needed is 16.7 ms. Some 
experimental results of using this processing 
sequence are shown in Chapter 9. 


8. Parallel architectures 

Yalamanchili et al. proposes a parallel archi- 
tecture in which the Hough space is divided into 
P partitions” and each processing element (PE) 


: N 
is required to increment Fa accumulators 


instead of Ng for each feature point. 

The parallel architectures in which each of 
the different @ coordinates is assigned a PE is 
an extreme case in which Ng equals P. We refer 
to this as completely parallel architecture. The 
completely parallel architectures that fall within 
the framework proposed by Yalamanchili et al. 
are the fastest, but they also require the most 
hardware resources. 

When the degree of parallelism increases, the 
percentage of circuitry dedicated to multiplica- 


Fig. 3—p calculation section. 


XY 
Edge extraction section FIFO ae 


PEn : Processor element 

TMCn: Table memory for x; cos 4@+n 
TMSn: Table memory for 4; sin 46°” 
HMzx : Histogram memory 


UO O UO O 
Output bus for 
histogram memory 


Peak detector 


6 address 
decoder 


Enable signals for HM, 


One Ng processor element is allocated to each @ coordinate. 
Neither multiplication nor transcendental calculation is 
required at run time. 


Fig. 4—Parallel architecture for Hough transform. 
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tion is the most important factor affecting 
economy. One set of parallel multipliers (12 
bits x 12 bits) consists of 1800 gates, so if 
Ne, is 512, a total of 920 K-gate are required 
for multiplication alone. Trigonometric func- 
tions, however, become less significant as the 
degree of parallel processing increases. This is 
because in completely parallel architectures, 
each PEn need hold only one set of cosmAg and 
sin nAg values. 

Here, we propose a completely parallel 
architecture that does not use any multipliers. 
Figure 4 is a block of diagram of this archi- 
tecture. 

Each x; and y; input can take 256 different 
values, so each x; cosnMg and Jy; sinnAg 
result for one PE can also have one of a possible 
256 values. This means that a 2 x 256 x 10 bit 
table can be used instead of the multipliers. 
Multiplication results were assumed to be 10-bit 
values after taking into account that the sum of 
both results would be quantized to nine bits. 

This architecture is advantageous from the 


point of view both of scale and of speed. First, if 


we assume a completely parallel architecture 
implemented with VLSI, we can compare the 
proposed system with a multiplier configuration 
from the standpoint of chip area. On one hand, 
a recent Fujitsu gate array design contains 
45000 gates in an 11.5 x 11.8 mm? area”, 
so if the above-mentioned multiplier of 12 x 12 
bits is configured into it, it will occupy 
5.43 mm”. On the other hand, Fujitsu’s 4 M-bit 
RAM with an access time of 70 ns measures 


21) If the above 


4.84 x 13.16 mm? in report 
table is implemented with the same technology, 
it will occupy 0.078 mm’. The area required by 
the lookup table is thus about 1/70 of that 
required by a multiplier with 1 800 gates. These 
values may contain some errors. To be com- 
pletely accurate, the multiplier should also have 
a circuit configuration suitable for 8-bit x 12-bit 
multiplication. The purpose here, however, is to 
estimate relative and qualitative trends of both 
systems. 

The operating speeds of systems based on 
table look-up architecture are faster because the 
table capacities are small and decoding time is 
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a) Hysteresis curve 


c) (p, 6) histogram 


d) Recognition results 


Fig. 5—Automatic measurements. 


much shorter than for 4 M-bit RAM chips. 


FUJITSU Sci. Tech. J., 24, 1, (March 1988) 


a) Magnetic disk head 


c) (p,@) histogram 
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d) Recognition results 


Fig. 6—Inspection applications. 


The proposed system has another advantage. 
It can be assembled using the most primitive 
type of commercially available microprocessors, 
such as those without multipliers. Each micro- 
processor references the table twice per feature 
point, performs one addition operation and 
increments the contents of the indicated 
memory address by the added value plus the dis- 
placement. The time required for these opera- 
tions can be in the order of a few tens of micro- 
seconds. 


9. Processing examples 
This chapter presents the experimental 
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results obtained using the model processor on 
three different images. Figures 5-a) to -d) and 
Figs. 6-a) to -d) show the results of using this 
processor in automatic measurement and inspec- 
tion applications. Figures 7-a) to -f) show the 
effects of using the modified complementary 
Hough transform. 

Figure 5-a) is a hysteresis curve of a thin- 
film head of a magnetic disk unit monitored by 
a curve tracer. Shown here is the upper half of 
the image plane. The main items to be measured 
are the upper and lower limits of the magnetic 
flux density (vertical axis) and the slope of the 
tangent line. The horizontal axis represents the 
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a) Mechanical part b) Picture based on extracted edges 


d) Modified complementary Hough transform 


e) Cells in Fig. 7-c) whose counts exceed one half the f) Cells in Fig. 7-d) whose counts exceed one half the 
peak height because of short line segments peak height because of short line segments 


Fi.g 7—Modified complementary Hough transforms. 


external magnetic field. Figure 6-a), taken head. Our objective was to detect any oblique 
through a microscope, shows the magnetic disk cracks or damage, which show up as the one 
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short and two long black lines in the figure. 

Figure 7-a) is an image of a mechanical part. 
Figures 5-b), 6-b), and 7-b) are based on extract- 
ed edge points. 

Figures 5-b), 6-b), and 7-b) show, respective- 
ly, the use of zero crossing after the Laplacian 
(4x4 Laplacian (2) in Table 1), enhanced 
vertical differentiation (4 x 4 enhancement (1), 
and the lower of 4x4 Sobel (1) in Table 1 
cascaded) using Sobel operation (4 x 4 Sobel (1) 
in Table 1). 

Figures 5-b), 6-b), and 7-b) respectively 
show the use of zero crossing after the Laplacian 
from the edge images. 

The recognition results were converted back 
to their original picture space coordinates and 
overlaid on the original pictures, as shown in 
Fig. 5-d) and 6-d). 

In these examples, thresholding and non- 
maximum suppression are used. The thresholds 
are 40 and 25, and the 5 x 5 and 7 x 7 windows 
are searched for maxima. 

Figure 7-d) shows the (p, @) space based on 
the modified complementary Hough transform. 
Applying the Laplacian after smoothing (4 x 4 
Laplacian (2) and 4 x 4 smoothing in Table 1) to 
Fig. 7-c) results in Fig. 7-d). 

Figures 7-c) to -f) provide evidence of how 
this method sharpens peaks and reduces sidelobe 
bias!®). Figures 7-e) and -f) show the cells whose 
counts exceed one half of the peak height as a 
result of the two parallel short line segments for 
Figs. 7-c) and -d), respectively. 


10. Conclusion 

We have developed dedicated straight-line 
detection hardware that uses the Hough trans- 
form method, and which is only slightly affected 
by noise and gaps even in complicated back- 
grounds. This hardware takes less than one 
second to produce straight-line parameters after 
video image capture. 

We have demonstrated the good recognition 
capabilities of the processor in automatic inspec- 
tion and measurement applications in factories 
and laboratories. 

We have proposed a highly parallel architec- 
ture that does not require any multiplication or 
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trigonometric calculations at run time. The 
advantages of this architecture have been dis- 
cussed from points of view of scale and speed. 

Future work will be directed towards in- 
corporating this processor into robot vision 
systems. 
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Reproducible Nb/AIO,/Nb 


Josephson Junctions 


@ Shiro Ohara @ Takeshi Imamura @ Shinya Hasuo 


(Manuscript received November 27, 1987) 


Nb/AIO,/Nb Josephson junctions were fabricated with different shapes and structures, and 
the influence of control lines, the size of the base electrode, contacts near the junctions, 
and junction size on junction characteristics was studies. 

Junctions fabricated by dc sputtering were found to be of high quality, regardless of shapes 


or structures. 


1. Introduction 

The Nb/AIO,/Nb Josephson junction has 
potential as an integrated Josephson circuit 
element because of its high quality and stability. 
The authors have developed the technology for 
fabricating Nb/AIO, /Nb junctions, and applied 
the junctions to logic circuits. High-speed opera- 
tion was confirmed for integrated Josephson 
circuits having several hundred gates. 

To utilize the Nb/AIO,./Nb junctions in cir- 
cuits, the junctions must be of high quality. This 
requires small subgap currents and sharp increas- 
ing currents at gap voltages. Critical current 
uniformity in a chip is also important. In study- 
ing the fabrication techniques for Nb/AIO, /Nb 
junctions the authors have clarified the ralation- 
ship between junction quality and deposition 
conditions”. In the course of this study, the 
authors encountered many difficulties in re- 
producing high-quality junctions. First, poor- 
quality junctions were frequently found on the 
same wafer with high-quality junctions. Second, 
the characteristics of high-quality junctions 
often deteriorated during subsequent processes, 
such as the formation of an insulation layer and 
Nb control lines. This phenomenon appears to 
be related to the size and shape of the junctions, 
base electrode, and control lines. This also seems 
to depend on the configuration of junctions and 
contacts. 
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Figure | compares the I-V characteristics 
observed in the same chip. Junction materials 
(Nb and Al) were deposited by rf magnetron 
sputtering. Figures l-a) and -b) show micro- 
photographs of junctions, and Figs. l-c) and 
-d) are their respective current-voltage charac- 
teristics. Figure l-a) is a pair connected in 
series consisting of a junction and a contact 
between the base electrode and the wiring layer. 
Figure 1|-b) is 6-um-diameter junctions connected 
in series. Despite the fact that they are on the 
same chip, their current-voltage characteristics 
are quite different. The quality parameter V,, is 
1.9 mV for the junctions in Fig. l-a) and 21 mV 
for those in Fig. 1-b). The critical current 
density J, for Fig. l-a) is five time that for 
Fig. 1-b). Here, Vm is defined as I, x Reg 
(/.: critical current, Rsg: subgap resistance at 
2mV). The only difference between these 
junctions is whether the contact is beside the 
measured junction. This suggests that junction 
characteristics are affected by the junction 
configuration. However, these deteriorated junc- 
tions could not be reproduced, and the causes 
of the deterioration are not clear. 

To study this irregularity in the junction 
characteristics, the authors fabricated Nb/ 
Al1O,/Nb junctions with different dimensions 
and configurations, and evaluated their charac- 
teristics. This time, junctions were deposited by 
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a) Junctions with contacts 


c) I-V characteristics for a); vert.: 2 mA/div, 
hor.: 0.5 mV/div. : 
(Vm = 1.9 mV, J.= 13000 A/cm*) 


——j 


b) Junctions without contacts 


d) I-V characteristics for b); vert.: 0.5 mA/div, 
hor.: 50 mV/div. 
(Vm = 21 mV, J, = 2700 A/cm?) 


Fig. 1— Microphotographs of junctions on the same chip and its I-V characteristics. 


dc magnetron sputtering, as this technique 
yields junctions of higher quality than those 
fabricated by rf magnetron sputtering!). (There- 
fore, dc magnetron sputtering is now used for 
circuit fabrication?’.) This method is also used 
because it is important to study the properties 
of junctions fabricated by methods currently 
used for integrated circuits. 
‘The authors examined: 
1) the influence of the control line above the 
junctions, 
2) the influence of contacts near the junctions, 
and 
3) the influence of the dimensions of the base 
electrode. 


2. Fabrication 

Figure 2 shows cross sections of the test 
devices. The junction structure consists of a Nb 
base electrode (B), an Al-AlO, tunneling barrier, 
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Fig. 2—Cross section of test deveices. 


and a Nb counter electrode (CT). These compo- 
nents were deposited on a wafer two inches in 
diameter in the same vacuum run. In addition 
to the junction electrodes, three metal layers 
were deposited: a Nb ground plane (GP), a Nb 
wiring layer (WR), and Nb control line (CL). 


FUJITSU Sci. Tech. J., 24, 1, (March 1988) 


They were insulated by SiO, except at the 
contact holes. All layers were deposited by 
sputtering. Nb and Al were deposited by dc 


Table 1. Thickness and etching gas 
Layer | Thickness (nm) Etching gass 
Gp 400 CF, + 20% O; 
I, 300 CHF3 + 15% O, 
B 250 | CF, + 10% 0, 
CT 200 | CF, +5% 0, 
I, 400 CHF; + 15% O, 
WR 600 | CF, +10% O, 
I; 800 | CFy +15% 0, 
CL 1 000 CF, + 5% O, 


me Se 


c) “Without CL” (Vm =41 mV, J, = 1400 A/cm?) 
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magnetron sputtering, and SiO, was deposited 
by rf magnetron sputtering. All layers were 
patterned by reactive ion etching (RIE). The 
surface of the B layer and the edge of the 
junction were oxidized by liquid anodization. 
The thickness and etching gas for each layer 
are listed in Table 1. 

Junctions having different dimensions and 
configurations were made on the same wafers. 
To study the influence of the I, and CL layers, 
the wafers were cut into two parts after the 
wiring layer (WR) was patterned. For measure- 
ment, one of them was diced into chips, and the 
other was diced after the subsequent deposition 
and patterning of the I, and CL layers. Chips 
diced before the deposition of the I, layer are 


b) “With CL” (Vj, =41 mV, J. = 1400 A/cm?) 


d) “With CL” (Vm =68 mV, J. =1490 A/cm?) 


5 


Fig. 3—Junction characteristics of one wafer and the other wafer (vert.: 0.2 mA/div, hor.: 1 mV/div). 
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termed chips ‘‘without CL” and chips diced 
after the CL layer was patterned are termed 
chips “with CL’’. 


3. Experimental results 
The junctions were evaluated by the param- 
eter V,,, and by critical current uniformity. 


3.1 Influence of the CL layer 

First, the influence of I, and CL layer 
deposition on junction characteristics was 
evaluated. Figure 3 compares the characteristics 
of junctions on chips “‘without CL” and those 
on chips “with CL’. The junction is 6.5 ym in 
diameter, and the control line is 2.5 um wide. 
Figures 3-a) and -b) show the change in the 
junction characteristics after the formation 
of the I, and CL layers in one wafer. Figures 3-c) 
and -d) show the change in the other wafer. 
Figures 3-a) and -c) represent the junction 
characteristics before deposition of I, layer. 
On one wafer, V,, decreases from 57 mV to 
41 mV after the formation of the I; and CL 
layers, and on the other wafer, V,, increases 
from 40 mV to 68mV. As a result, we could 
not observe any relationship between the quality 
of the junction and the deposition of I, and CL 
on the junctions. Variation in V,, shown in 


Table 2. Current density 


(A/cm?) 


Control line 
width (um) 


j ; 4.0 
unction 
diameter (um) 


3.5 ias0 | 1260 | -— 
4.5 1510 | 1380 | 1380 
6.5 1450 | 1390 | — 


Table 3. Quality parameter V,, 


(mV) 


Control line 
width (tm) | 9 25 | 4.0 
Junction 
diameter (um) 
3.5 63 52 _ 
4.5 44 42 42 
6:5 62 53 - 
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Fig. 3 should be considered due to Vy, scatter- 
ing on the wafer. For the critical current density, 
J~, changes induced by the formation of I, and 
CL layers were also small. This was the same 
as the intrinsic J, spread on a wafer. 

The critical current density and V,, measured 
for different junction areas and control line 
widths are listed in Tables2 and 3. Three 
different junction areas were prepared ranging 
from 3.5 um to 6.5 wm in diameter. The control 
lines were 2.5 um and 4.0 ym wide. A junction 
where the deposited CL was completely removed 
was also prepared on the chip. This is referred to 
as having a CL width of Oywm. No distinct 
relationship between the junction characteristics 
and control line widths and/or junction diam- 
eters were observed. 

The influence of CL layers on the J, uni- 
formity was also evaluated for 100 junctions 
connected in series. The results are shown in 
Table 4. The diameter of the junction area is 
6.5 um-diameter. The control line is O um, 
2.5 um, and 22.5 um wide for 25-junction 
chains. The J, spread ranges from 5.6 percent 
to 7.4 percent. This does not indicate dependence 
on the width of the control line. It is almost the 


Table 4. Control line width and critical current spread 
(junction dia.: 6 um) 


Control line width I, spread Number of 
(um) (%) junctions 

0 6.0 25 

2.5 7.4 25 

22.5 5.6 25 

= 6.4 100 


Table 5. B layer sizea and junction characteristics 


B layer size Current density Vm 

(um? ) (A/cm? ) (mV) 
3 250 650 ae 
1750 680 51 
1350 740 49 
1000 700 49 
700 680 EW 
450 690 $1 
250 730 52 
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same as the value of 6.4 percent measured for 
the 100-junction chain without the CL layer. 
It is reasonable, therefore, to conclude that the 
formation of the I; and CL layers has no sig- 
nificant influence on J, spread. 


3.2 Influence of base electrode size 

We evaluated the characteristics for junctions 
with different base electrode sizes. Other junction 
dimensions are the same for all samples. The 
critical current density and V,, are listed in 
Table 5. The area of the base electrode varies 
from 250 um? to 3250 um?, and the junction 
area is 4.5 um in diameter. The current density 


a) With contacts 


30 um 
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and V,, are quite constant for all junctions, 
and no relationship with the base electrode 
size was observed. 


3.3 Influence of contacts 

When the contact between the base electrode 
and the wiring is placed near a junction the 
characteristics of the junction often deteriorated 
in the same way as shown in Fig. 1. To study the 
influence of the contact, pairs of junctions and 
contacts were prepared in series, and compared 
with the simple junctions connected in series. 

Figure 4 shows microphotographs of the 
test devices. The junction is 6.5 um in diameter, 


b) Without contacts 


Fig. 4—Junction array for evaluation the influence of contacts. 


a) With contacts (V;, = 40 mA, J, spread: 11 percent) 


b) Without contact (Vm = 40 mV, /, spread:7.7 percent) 


Fig. 5—I-V characteristics of 100 junctions connected in two series (vert.: 0.1 mA/div, hor.: 50 mV/div). 
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and the contact is 7 um square. Figure 5 shows 
I-V characteristics for 100 junctions connected 
in series, Figs 5-a) with contacts and -b) without 
contacts. The value of Vj was evaluated at 
40 mV for both samples. Here, an /, distribution 
median was used to calculate V;,. The /, spread 
was slightly different, 11 percent for Figs. 5-a) 
and 7.7 percent for -b). In Fig. 4a), a few 
junctions with very large 7, were observed in the 
I-V characteristics, and they apparently enlarged 
the J, spread. If they are omitted, /, spread is 
thought to be almost constant. Thus, the authors 
concluded that the influence of the contact near 
the junction is not observed in the junction 
characteristics. 


4. Conclusion 

Junctions with different sizes and configura- 
tions were fabricated and evaluated to determine 
the influence of the control line, the base 
electrode, and the contact on junction charac- 
teristics. No distinct differences in the junction 
quality and the J, spread were observed among 
junctions despite the different shapes and 
structures. 

The inconsistencies in the production of 
high-quality junctions, as shown in Fig. 1, were 
observed frequently during the development 
of fabrication techniques. When junction elec- 
trodes were made by rf sputtering and the 
maximum V,, was no more than 25 mV, it was 
observed that the junction characteristics seemed 
to depend on the junction shapes. This time, 
junctions were fabricated by dc _ sputtering 
because this method yields superior junction 
characteristics. The quality parameter, V,, 
consistently exceeded 40 mV in all shapes and 
constructions tested. No influence of the shape 
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or structure on the junction characteristics 
could be observed. Junctions fabricated by 
dc sputtering always have high-quality charac- 
teristics, regardless of their shapes and structures. 
The new fabrication technology is suitable for 
Josephson integrated circuits consisting of 
different shapes and structures. Deterioration 
in the junction characteristics seems to be 
prevented by increasing the value of V,,. How- 
ever, the question of why the deterioration 
occurred and why dc sputtering suppresses this 
remain unanswered. Considering that the deterio- 
ration could be prevented by changing the 
method of film deposition, the internal film 
stresses or the conditions around the junction 
interface seem to cause the deterioration. Why 
this is the cause is not yet clear. The determina- 
tion of the cause of the deterioration in charac- 
teristics is important to fabricating stable and 
reproducible high-quality Josephson junctions. 
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Resonant-Tunneling Hot Electron 
Transistors (RHETs) Using 
an InGaAs/In(AlGa)As Heterostructure 
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This paper describes the RHETs using an InGaAs/In (AlGa) As heterostructure. A common- 
emitter current gain of 10.0 was measured at 77 K for a 100 nm base RHET. The measured 
current gain of 25 nm base RHETs reaches as high as 25.0, which is about five times larger 


than that of GaAs/AlGaAs RHETs. 


The collector current peak-to-valley ratio reached 15.0 for the same device. 


1. Introduction 

In 1985, we proposed and fabricated three- 
terminal resonant-tunneling devices, named the 
RHET  (Resonant-tunneling Hot Electron 
Transistor), which used a GaAs/AlGaAs hetero- 
structure!), We have demonstrated various 
attractive features of RHETs for memory and 
logic applications??. However, the common- 
emitter current gain in GaAs/AlGaAs RHETs 
was as low as 5 for a 25 nm base, which is not 
high enough for practical use. 

The relatively low current gain in GaAs/ 
AlGaAs RHETs is believed to be due to the 
narrow I-L valley separation energy of GaAs 
(0.30 eV). On the other hand, in InP and 
Ino.53Gag.47As, the I-L_ valley separation 
energy is considered to be about 0.55 eV. 

In InAs, the I’-L valley separation energy is 
considered to be about 0.78 eV. So, these 
materials should work well as the base layer for 
RHETs. Here, we tried to fabricate InGaAs base 
RHETs, which can be easily grown on InP 
substrate using MBE without lattice mismatch. 


2. Experiment 

Figure 1 is a schematic cross section of the 
InGaAs/In(AlGa)As RHET we _ fabricated. 
InGaAs and In(AlGa)As layers were grown on 
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the InP semi-insulating substrate using pulsed 
molecular beam methods*). The device uses 
a quantum well which acts as the hot electron 
injector. This 3.81 nm 
Ing.53Gap.47As layer sandwiched between two 
layers of 4.4 nm Inp.53 Alo.47As. We fabricated 
two types of RHETs with different base widths. 
The thicknesses of the base layers were 25 nm 
and 100 nm. The carrier concentration of the 
base layer was 1 x 10!§ cmv?. The thickness of 
the Inp.52(Alo.5 Gag.5)o.4gAs collector barrier 
was 200 nm with a 6 nm compositional stepwise 
graded layer to reduce the quantum mechanical 
reflection at the collector barrier. 

To investigate the effects of the graded 
layer, we also fabricated 25 nm base RHETs 
without the graded layer. 

All the current-voltage characteristics were 
measured at 77 K. 


consists of a 


3. Results 

Figure 2. shows the common-emitter col- 
lector current-voltage characteristics measured at 
77K for a 100 nm base RHET. A common- 
emitter current gain of 10.0 was obtained at a 
collector-emitter voltage of 2V and a _ base 
current of | mA. The current gain was about 15 
times larger than that of the 100nm_ base 
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n-InGaAs 


InAlAs (4.4 nm) 
InGaAs (3.81 nm) 
InAlAs (4.4 nm) 


i-InAlGaAs 
(200 nm) 


InP substrate 


Fig. 1—Schematic cross section of InGaAs/In(AlGa)As 
RHET. 


100 nm base 


: 200 wA/step 


Collector current (mA) 


0 0.5 1.0 15 2.0 2.5 


Collector-emitter voltage (V) 


Fig. 2-Common-emitter collector current-voltage 
characteristics for the 100 nm-base RHET. 


GaAs/AlGaAs RHETs®. Using thermionic 
emission current measurement, the conduction 
band discontinuity at the collector barrier was 
found to be 0.28 eV, which is about half that of 
the InAlAs ternary barrier (0.55 eV). The 
increased current gain of InGaAs/In(AlGa)As 
RHETs compared to that of GaAs/AlGaAs 
RHETs is believed to be due to the larger T-L 
valley separation energy of InGaAs compared to 
that of GaAs. 

Figure3 shows the common-emitter 
collector current-voltage characteristics for the 
25 nm-base RHET with a graded layer. 
A common-emitter current gain of 18 was 
obtained at a collector-emitter voltage of 2 V, 
which is about twice that of the 100 nm-base 
RHET. Also a high current gain of 25 was 
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25 nm base 


Ig: 100 nA/step 


WA 


eee 


Collector-emitter voltage (V) 


Collector current (mA) 


Fig. 3—Common-emitter collector current-voltage 
characteristics for the 25 nm-base RHET. 


10 


Jp = 2.8 X 104 A/cm? 
Je/Jy = 15.0 


25 nm base 


Collector current (mA) 


0 
0 0.5 1.0 1.5 2.0 


Base-emitter voltage (V) 


Fig. 4—Collector current as a function of base-emitter 
voltage for the 25 nm-base RHET. 


obtained at a collector-emitter voltage of 3 V, 
which is about five times larger than that of 
GaAs/AlGaAs RHETs. Figure 4 shows the 
collector current as a function of the base- 
emitter voltage for the same device. The 
collector current exhibits a pronounced peak 
due to resonant-tunneling at a_ base-emitter 
voltage of 1V. The peak collector current 
density is 2.8 x 10* A/cm? and the peak-to- 
valley ratio reached 15.0. The increased current 
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Non-graded 


W, = 25nm 


: 200 nA/step 


Collector current (mA) 


0 0.5 1.0 1.5 2.0 2.5 3.0 


Collector-emitter voltage (V) 


Fig. S-Common-emitter collector current-voltage 
characteristics for the 25 nm-base RHET 
without graded layer. 


density and peak-to-valley ratio is considered to 
be due to the large tunneling probability and 
high barrier height of the emitter InAIAs barrier 
compared to that of the AlGaAs barrier. The 
high current gain and peak-to-valley ratio should 
contribute to the stable operation of the 
circuits. The high peak current density should 
contribute to the high-speed operation of the 
RHET because it reduces the charging time for 
the emitter and collector capacitances. 

Next, the effect on the current gain of the 
compositional graded layer in the collector 
barrier was investigated. We fabricated two kind 
of 25 nm-base RHETs. One has a 6 nm three 
stepwise compositional graded layer at the base 
and collector barrier interface (with a graded 
layer). The other has no graded layer (without 
a graded layer). Figure 5 shows the common- 
emitter collector current-voltage characteristics 
for the 25 nm-base RHETs without a graded 
layer. A common-emitter current gain of eight 
was obtained at a collector-emitter voltage of 
2 V and a gain of ten was obtained at a col- 
lector-emitter voltage of 2.5 V, which is smaller 
than that of the devices with a graded layer. 

Figure 6 shows the common-emitter current 
gain as a function of collector-emitter voltage 
for the graded and non-graded RHETs. Current 
gain was measured at the base current range 
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Graded g 
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oe 
Z Non-graded 
—— 


Current gain 


Vee (V) 


Fig. 6—Common-emitter current gain as a function of 
collector-emitter voltage for graded and 
non-graded RHET. 


Non-graded 


Barrier height (eV) 


Voltage (V) 


Fig. 7—Collector barrier height as a function of applied 
voltage for non-graded RHET. 


from 200 uA to 400 uA. The current gain for 
the graded RHET is lower than that of the non- 
graded RHET at a lower Vcg than 1.5 V (Vcp 
of about 0.9 V). 

Figure 7 shows the barrier height as a func- 
tion of the applied voltage for non-graded 
collector barrier RHETs. The measured results 
are shown for the forward and reverse bias. The 
barrier height decreased as the applied voltage 
was increased. This figure clearly shows that 
there is a difference between the forward and 
reverse bias barrier heights and that the 
difference was almost constant between 0.01 V 
and 1 V. This difference can be explained as the 
difference of Ec-Ep (33 meV) due to the 
difference in the doping concentration for the 
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Fig. 8—Collector barrier height as a function of applied 
voltage for graded RHET. 


base layer and collector layer (Np =1 x 
10!8 cm™3 for the base layer, Np = 5 x 10!7 cm’? 
for the collector layer). Figure 8 shows the 
barrier height as a function of the applied 
voltage for the 25 nm base RHET with the 
graded layer. The difference in the barrier 
heights between the forward bias and the reverse 
bias is increased as the applied voltage is in- 
creased. This increase in the barrier height 
difference is believed to be caused by the low- 
ering of the graded barrier with the applied 
voltage. Comparing the results for non-graded 
barrier RHETs, the barrier height lowering with 
an applied voltage of 1 V is calculated to be 
about 30meV, which agrees well with the 
theoretically calculated barrier lowering. The 
initial barrier height lowering should contribute 
to the difference in the breakdown voltage in 
the common-emitter configuration. The charac- 
teristics shown in Fig. 6 can be well explained as 
the increase in current gain caused by the barrier 
height lowering at positive collector-base voltage 
for the graded RHETs. The effect of quantum 
mechanical reflections on the current gain was 
not clearly determined in this investigation. 

Next, we estimated the device transit time 
for the InGaAs/In(AlGa)As RHETs. The total 
RHET delay time can be represented using 
Equation (1). 


Tec = TE + TR TTC * TEC- 


Here, te is the emitter depletion-layer 
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charging time, 7g is the base transit time, Tc is 
the collector depletion layer transit time, and 
Tcc is the collector charging time. It is interest- 
ing to compare the total delay time of RHETs 
with that of GaAs/AlGaAs HBT (Heterojunction 
Bipolar Transistors). The same device size was 
assumed for RHETs and HBTs. The emitter/base 
junction was 2 um x 5 wm and the base/collector 
junction was 4 um x 7 um. Cheng et al. reported 
a high cut-off frequency of 67 GHz. The 
individual delay time was 0.6 ps for Tp, 0.6 ps 
for Tg (80 nm base), 0.7 ps for tc, and 0.5 ps 
for tcc. The total delay time was 2.4 ps”. 

For the InGaAs/In(AlGa)As RHETs, the 
base width was assumed to be 100 nm as was the 
thickness of the collector barrier. 

The base transit time was calculated assum- 
ing ballistic transport in the base layer. The 
injection energy was supposed to be 0.6 eV, so 
the electron velocity was calculated to be about 
1.24 x 108 cm/s. The base transit time was 
calculated to be about 0.08 ps”. 

The collector depletion layer transit time 
was calculated to be about 0.01 ps, considering 
the velocity change at the base/collector barrier 
interface. The collector charging time was 
calculated to be about 0.06 ps. Then the total 
delay time, excluding the emitter depletion layer 
charging time, was calculated to be about 
0.24 ps, which is about 1/7 less than that of the 
HBT. These results indicate the advantage of hot 
electron transistors for collecting only ballistic 
or near ballistic electrons. 

The emitter delay time is related to the built 
in time, ¢;, to achieve the resonant-tunneling 
state and the time, f2, for electrons passing 
through the resonant-tunneling barrier. The time 
ty was calculated using h/AE). AEo was cal- 
culated to be 0.51 meV for the barrier used in 
this paper, so f2 was calculated to be 1.3 ps. The 
time f; was considered to be the charging time 
for an ordinary parallel plate. The thickness of 
the capacitance was considered to be roughly 
the sum of the thicknesses of the barrier and the 
depletion layer. The conductance was calculated 
to be about the peak current divided by half of 
the peak voltage, which agrees well with the 
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measured results. The calculated value of tf; was 
2.4 ps. The simple sum of ¢; and f was as large 
as 3.7 ps. Actually, the total emitter delay time 
was considered to be smaller than f; + f2. If we 
assume that the total emitter delay time is the 
larger of ¢; and f2, the emitter delay time can be 
calculated to be about 2.4 ps. In this case, the 
total delay time was calculated to be 2.64 ps 
which is almost equal to that of HBTs. 

If the thickness of the emitter barrier was 
decreased, Eg increased and fz decreased. A fy 
also decreased because of the increased current 
density and therefore increased conductance. 
The decreased barrier thickness (3.52 nm) and 
barrier width (3.22 nm) cause a peak current 
density of 4 x 10° A/cm? and the delay time f; 
will be about 0.62 ps (t7 = 0.21 ps). In this case, 
the total delay time was calculated to be about 
1 ps and the high cut-off frequency of 160 GHz 
will be achieved. These calculations indicate that 
InGaAs/In(AlGa)As RHETs have great potential 
as high-speed devices. 


4. Conclusion 

We have fabricated InGaAs/In(AlGa)As 
RHETs. A common-emitter current gain of 10 
was obtained for the 100nm_ base RHETs, 
which is about 15 times larger than that of 
GaAs/AlGaAs RHETs. A high current gain of 25 
was obtained for the 25 nm base RHETs with 
a compositional stepwise graded layer. The 
collector current peak-to-valley ratio reached 15. 
The increased current gain compared to that of 
GaAs/AlGaAs RHETs is considered to be caused 
by the large I-L valley separation energy of 
InGaAs compared to that of GaAs. Compared to 
the current gain for the 25 nm base RHETs 
without a graded layer, the barrier height lower- 
ing at the applied voltage contributed to the 
increase in current gain. 

The device delay time was estimated for 
RHETs compared to that of GaAs/AlGaAs 
HBTs. The estimated delay time was about | ps 
with a cut-off frequency of 160 GHz. We found 
that the RHETs have great potential as high- 
speed devices. 
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The obtained high current gain and high 
peak-to-valley ratio indicated that the InGaAs/ 
In(AlGa)As RHETs are very promising for high 
speed and stable operation of future integrated 
circuits. 
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A new type of crystal phase change was found in selenium-indium-antimony alloy film, 


which is a unique material developed by Fujitsu. The phase change occurs reversibly between 


two different crystal states when the film is irradiated with laser pulses, resulting in different 


optical reflectivities. An erasable optical disk using this film has good carrier-to-noise-ratio of 


45 dB, moderate recording sensitivity and long media life, but erasable cycles still vary from 


between 10° and 10° times at present. Though there remain some problems to be solved, 


the disk is nearly ready for practical use, and as such is very promising. A proposed model 


of the phase-change mechanism is also discussed. 


1. Introduction 

Optical disks have a high potential for large- 
volume information storage. Using a diffraction- 
limited focused laser beam (1 wm), it easily 
enables a storage density of 107-10® bit/cm?’, 
which is ten times larger than the storage density 
of the most advanced hard magnetic disks. 
Different from magnetic disks, optical disks 
store information on the thin film which is 
sandwiched between two. relatively thick 
(1.2mm) substrates. Because writing and 
reading are performed through the transparent 
substrate, dust on the substrate surface causes 
very few disk errors. Thus, optical disks can 
serve as easy-to-handle and high-performance 
file memories for a large computer system as 
well as for a small computer used in an office 
or at home. 

Read-only type disks, from which only 
prerecorded information can be read, are already 
available in the form of video disks or CDs. 
As for write-once disks, on which information 
can be recorded by users, Fujitsu has already 
been shipping write-once disk drives (FACOM 
6441 and 6442). Fujitsu. has also been 
developing erasable optical disks, focusing on 
two types of erasable disks: magneto-optical 
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disks and phase-change disks. This paper is 
about phase-change erasable disks. 

The key to developing a _ phase-change 
erasable disk is to find a new material. The 
authors found out that a _ selenium-indium- 
antimony (SeInSb) alloy film shows reversible 
phase change when irradiated with laser light, 
and is thus promising as a recording material 
for erasable disks’. The authors also have 
found that an alloy film of indium and 
antimony (InSb) has the same erasability”. 
This paper describes the results of experiments 
evaluating the characteristics of the SeInSb disk 
and analyzing light-irradiated micro-bits of the 
alloy films. The paper also describes a proposed 
model of the phase-change mechanism of the 
alloy films®~®, 


2. Principle of writing/erasing method on phase- 

change media 

The writing procedure on _ phase-change 
media changes the crystal state of the thin film 
by irradiating it with a laser beam. Phase 
transition is caused by the thermal effect of the 
irradiated laser light. No other external energy, 
such as an external magnetic field is used. 
Generally, two different irradiating conditions 
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Overcoating layer 
(Photopolymer : 30 um) 


Recording layer 
(SeInSb: 100 nm) 


Under layer 
(ZnS: 20 nm) 


Groove-replicating layer 
(Photopolymer : 50 wm) 


Substrate 
(Glass: 1.2 mm) 


Laser beam 


Fig. 1—Cross section of SeInSb medium. 


are used. One is to irradiate at high power for 
a short duration. The film changes to a liquid 
state in a short time and then solidifies quickly. 
The cooling rate is so high that the atoms cannot 
obtain energetically stable positions. The other 
condition is to irradiate at lower power for 
longer duration. In this case, cooling rate of 
melted film is so low that the atoms have time 
to recover energetically stable positions. 

Tellurium alloy films have been developed 
as erasable phase-change materials» ®. In these 
films, information was recorded using a phase 
change between an amorphous state and a 
crystal state of the films. The newly developed 
alloy films, SeInSb and InSb films, have a 
different phase-change mechanism. They use 
a phase change between two different crystal 
states. Information is recorded via a transition 
between these two stable crystal states. Laser- 
irradiating conditions are nearly the same as 
with the tellurium alloy films. 


3. Experiment 
3.1 Fabrication of alloy films 

SeInSb and InSb alloy films were deposited 
on a substrate by a vacuum. evaporation 
technique. Vapor sources of selenium, indium 
and antimony were heated independently. 
The three vapors were mixed together on 
rotating substrate and formed a thin alloy 
film. The evaporation rate of each vapor was 
controlled independently, using a  quartz- 
crystal oscillator to strictly control the film 
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Fig. 2—Photograph of SeInSb optical disk. 


composition. 

A cross section of the medium is shown 
in Fig. 1. The substrate is glass which is coated 
with UV polymerized photopolymer (2P)”. 
Guiding grooves for the tracking  servo- 
mechanism were formed on the surface of the 
2P layer. Between the 2P layer and the recording 
layer (alloy film), zinc sulfide was deposited to 
reduce the reflectivity of the medium and this 
increased writing sensitivity. The composition 
of the alloy film was 20 at% selenium, 35 at% 
indium and 45 at% antimony for SeInSb, and 
40 at% indium and 60 at% antimony for InSb. 
Finally, the top surface was coated with another 
2P layer. The purpose of the overcoating was 
not to protect the alloy films from oxidation, 
but rather to prevent hole formation in the alloy 
film when it is irradiated with laser light. 
Figure 2 is a photograph of a 130-mm-diameter 
SeInSb optical disk. 


3.2 Evaluation of write/erase characteristics 

3.2.1 Static recording 

SeInSb medium was irradiated with a 
focused light beam from a laser diode. In this 
case, the medium was not rotated. The diameter 
of the beam was 1 wm. The wavelength was 
830 nm. The position of the objective lense 
which focused the beam was automatically 
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Fig. 3—Reflectivity change by cyclic irradiation. 


controlled to focus the beam on the film in 
such a way that the beam diameter reached 
its minimum at the surface of the alloy film. 

Between irradiations for writing and erasing, 
a weak continuous power was applied and 
reflected light from the small area of the film 
was detected, which gave the reflectivity of the 
small irradiated area of the medium. As the laser 
light was irradiated repeatedly on the medium, 
the reflectivity changed as shown in Fig. 3. 
At first, as-deposited alloy film had low reflec- 
tivity. Then, after an erasing beam was applied, 
the medium reached a new state of slightly 
higher reflectivity. The next irradiation with 
a writing beam moved the medium to a still- 
higher state of reflectivity. The next irradiation 
with an erasing beam returned the medium to 
the low-reflectivity state. After that, cyclic 
writing and erasing changed the reflectivity 
reversibly between the high and low states, 
respectively. The reversible change of reflectivity 
was confirmed beyond 10° cycles. 

Figure 4 shows irradiating conditions 
necessary to make transitions between the two 
crystal states. An irradiation at high power for 
a short duration causes the transition from 
crystal state II to I. On the other hand, an 
irradiation at low power for a longer duration 
causes the reverse transition from [| to Il. 

3.2.2 Dynamic recording 

The irradiating method for a rotating disk 
is somewhat different from that for a stationary 
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Fig. 4—Irradiating conditions for changing reflectivity. 


disk. Since each point of a rotating disk is 
irradiated only for the time during which the 
point traverses the laser beam, the rotation 
speed and the beam diameter determine the 
duration of irradiation. As was described in 
subsection 3.2.1, the duration of an erasing 
pulse is longer than that of a writing pulse. 
So, in this experiment, in order to obtain 
two different pulse durations at a fixed rotation 
speed, a two-beam optical head was used. 
One of the beams was a conventional circular- 
shaped beam for writing and reading. The 
other beam was an elliptical-shaped beam 
for erasing. The long axis of the erasing beam 
was aligned in the direction of the data tracks. 
The diameter of the circular-shaped beam 
was 0.9 um. The elliptical-shaped beam was 
3 wm long and | wm in width (i.e. perpendicular 
to the track direction). 

The rotation speed was fixed at 900 rpm. 
The linear velocity was about 6 m/s. A signal 
of 3 MHz was recorded. The necessary peak 
power of the modulated laser pulse was 8 mW 
on the medium surface. The bit length was 
1 wm. The same signal was erased by irradiating 
an elliptical-sshaped DC laser beam on the 
recorded track. The necessary power was 
12 mW. Overwriting was also performed. The 
prerecorded signal was erased by the elliptical 
beam and just after that a new signal was 
recorded on the same track during the same 
revolution of the disk. 
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Recorded tracks 
(six tracks) 


Erased tracks 
(two tracks) 


10 wm 


Fig. 5—Photomicrograph of recorded and erased tracks. 


Considering the linear velocity of 6 m/s, 
laser power and pulse duration are 8 mW, 170 ns 
for writing and 3 mW, 500 ns for erasing. These 
conditions coincide well with the data for static 
recording in Fig. 4. 

A photomicrograph of the recorded tracks 
and erased tracks is shown in Fig.5. Bits of 
1 wm can be seen clearly. The carrier-to-noise 
ratio was measured to be 45 dB. The measuring 
bandwidth was 30kHz. After erasure, the 
3 MHz prerecorded signal was reduced to 
—40 dB of the initial level. Figure 6 shows the 
waveforms of the read-out signals of the 
recorded track and the erased track. 

An experiment of cyclic writing and erasing 
was performed. It was observed that the carrier- 
to-noise ratio reduced gradually as number of 
cycles increased. The maximum number of 
erasable cycles varied between 10° and 10° 
cycles from disk to disk. As was determined by 
static recording, the reversible change of reflec- 
tivily exceeded 10° times. In this experiment, 
the amplitude of the read-out signal, which 
corresponds to the reflectivity change in the 
static experiment, did not change as the number 
of cycles increased. But the noise, which could 
not be observed in the static experiment, 
increased with the number of cycles. 


3.3 Analysis of bits 

3.3.1 Static recording 

The irradiated SeInSb alloy film was peeled 
off from the substrate and observed by trans- 
mission electron microscope. The bright field 
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a) Recorded track 


500 ns/div 


b) Erased track 


Fig. 6 —Wave forms of read-out signals. 


images of the irradiated area (bits) and diffrac- 
tion patterns are shown in Fig. 7. The outside 
of the irradiated area, which remained in the as- 
deposited state, was amorphous. From the 
diffraction patterns it was clear that the ir- 
radiated bits were an assembly of small crystals. 
The difference between a written bit and an 
erased bit was not so definite, but the grain size 
of the central part of a written bit seemed to be 
bigger than that of an erased bit. Accordingly, 
the reflectivity change shown in Fig. 3 can be 
explained in terms of the crystalline states; 
the initial low-reflectivity state corresponds 
to the amorphous state. Laser light irradiation 
changed the film from the amorphous state to 
the crystal state. The film had two kinds of 
stable crystal states with different reflectivities. 
The cyclic irradiations with writing and erasing 
pulses caused reversible transitions between 
these crystal states, resulting in reversible 
reflectivity changes. 

The InSb alloy medium showed the same 
reflectivity change and had same reversible 
cycles as the SeInSb medium shown in Fig. 3. 
Figures 8-a) and 9-a) are transmission electron 
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a) Written bit 


b) Erase bit 0.5 ym 


Fig. 7—Transmission electron micrographs of SeInSb 
film. 


micrographs of InSb alloy film. In this case, 
the alloy film was crystallized during vacuum 
evaporation, so the as-deposited initial state 
is crystalline. It must be noted that the 
diameters of irradiated areas in Figs. 8-a) and 
9-a) are smaller than in Fig. 7. The reason is 
thought to be that in an amorphous film, heat 
diffusion from the irradiated area to the outer 
amorphous area results in an excess crystallized 
area around the bit. The crystallization during 
deposition can be also applied to SeInSb film, 
which is effective to obtain small bits. 

Compared with SelInSb, InSb is simpler. 
For this reason, the authors conducted a 


64 


b) Results of dark field 
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Fig. 8—Transmission electron micrographs of the written 
bit on InSb film. Identified crystals are marked 
b). 


thorough investigation of InSb film in order 
to clarify the mechanism of the crystal phase 
change. From the analysis of the diffraction 
patterns, the outside of the irradiated bit of 
the InSb film was shown to be a mixture of 
very small Sb crystals and small InSb compound 
crystals. In this film a clear difference exists 
between written bits and erased bits. In the 
written bit { Fig. 8-a)}, the center part is 
composed of large crystals. In the erased bit 
| Fig. 9-a)} , however, the center part is an 
assembly of small crystals. At the outer part 
of the both bits, crystals of intermediate size 
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Fig. 9—Transmission electron micrographs of the erased 
bit on InSb film. Identified crystals are marked 
b). 
Small crystals of Sb and InSb in the central part 
are not marked deliberately because of 
complications. 


are observed. From the analysis of the diffrac- 
tion patterns, it was concluded that all the 
crystals in both bits were composed of cubic 
InsgSbs9 intermetallic compound (zinc blend 
structure) and rhombohedral Sb. No other 
crystal structure was found. 

Dark field observation was also performed. 
Images of Sb crystals were formed by electrons 
that were diffracted from one of Sb crystal 
planes. In the same manner, images of InSb 
compound crystals were formed. Figure 8-b) is 
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the combined result of the dark field observa- 
tions. It suggests that the large crystals at the 
center part were Sb crystals, and the crystals 
of intermediate size at the outer part of the bit 
were InSb compound crystals. 

As for the erased bit, there were no large 
Sb crystals, but Sb and InSb compound crystals 
were distributed uniformly in the center part of 
the bit | Fig.9-a)}, which was determined 
from the analysis of the diffraction patterns. 
Result of the dark field observations is shown 
in Fig. 9-b), in which only InSb crystals at the 
outer part of the bit are marked. (Small crystals 
of Sb and InSb compound in the central part 
are not illustrated because of complications.) 

From the Figs. 8 and 9, it was concluded 
that the reversible change in InSb alloy film 
is related to the distribution of Sb crystals 
and InSb compound crystals. In SeInSb alloy 
film, too, the origin of the reversible change is 
thought to be change in distribution of Sb 
crystals and InSb compound crystals. The reflec- 
tivity change can also be explained in terms of 
the distribution change. It is known that reflec- 
tivity of bulk Sb and InSb compound is 
70 percent and 40 percent, respectively. So, 
the written bit, which has large Sb crystals at 
the central part of the bit, shows high reflectivi- 
ty. On the other hand, the erased bit, in which 
the central part is relatively poor in Sb crystals, 
has low reflectivity. 

3.3.2 Dynamic recording 

The recording layer was peeled off from the 
substrate and observed by transmission electron 
microscope. Figure 10-a) is the transmission 
electron micrograph of a bit array recorded 
on a rotating disk. This alloy film was also 
crystallized during vacuum deposition, because 
in the crystallized film smaller bits can be 
recorded than in amorphous film, as described 
in Subsec. 3.3.1. The recorded bits were located 
along the groove. At the center of each bit, large 
crystals appeared, which are assumed to be Sb 
crystals. 

Figure 10-b) is the transmission electron 
micrograph of an erased track. The bit patterns 
in Fig. 10-a) have disappeared completely, 
but along both sides of the track erased by 
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a) Recorded tracks 


b) Erased track 


Fig. 10—Transmission electron micrograph of recorded tracks and erased track on a rotating disk. 


beam-irradiation, crystals of intermediate size 
were present, and are assumed to be InSb 
compound crystals. These photographs are 
basically the same as the transmission electron 
micrographs of static recording bits shown 
in Figs. 8 and 9. 


3.4 Estimation of media life 

SeInSb and InSb alloy films are chemically 
stable. The oxidization of both films has been 
monitored by measuring the optical trans- 
mission. It has been confirmed that optical 
transmittances of both films do not change 
for more than four months under the 
atomospheric condition of 60 °C and 90 percent 
relative humidity” 

In this experiment, stability of the crystal 
phases was examined. That is, the storage time 
of recorded information was examined by an 
accelerated aging test. If the written state 
(crystal I) is energetically unstable, information 
will be lost during an aging test. The media 
life was evaluated in terms of carrier-to-noise 
ratio. A signal of 3 MHz was recorded on both 
SeInSb and InSb disks and their carrier-to- 
noise ratio were measured after certain intervals 
under the testing atmosphere. The humidity 
was fixed at 90 percent relative humidity, 
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Fig. 11—Arrhenius plot of media life vs. aging 
temperature. Life was defined end at the time 
that C/N was reduced by 3 dB from the initial 
value. 
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Fig. 12—Equilibrium phase diagram of InSb alloy system. 
Arrows show cooling processes of the written 
bit and erased bit. 


and only temperature was varied (60-120 °C). 

Life of the disk under the measured condi- 
tions was defined to be the time interval 
between the initial time and the specific time at 
which the carrier-to-noise ratio dropped to 
—3 dB from the initial value. The results are 
shown in Fig. 11. Extrapolation from the high- 
temperature region to 25 °C suggests that the 
life of an InSb disk is about 20 years at 25 °C 
and 90 percent relative humidity. The life of 
an SeInSb disk is estimated to be more than 
50 years. Selenium itself does not have a direct 
influence on the mechanism of the reversible 
phase change, but it has an important role in 
increasing chemical stability and data storage 
time. 


4. Discussion 
4.1 Model of phase-change mechanism 

Based on the observations described in 
section 3.3 and using equilibrium phase diagram, 
the authors have proposed a model of the 
writing and erasing mechanisms in SeInSb and 
InSb alloy films. Figure 12 is the phase 
diagram of the In-Sb binary alloy system. The 
authors have found that the alloy film in which 
Sb concentration is between 50 at% and 69 at% 
shows a_ reversible reflectivity change when 
irradiated with a laser pulse. This is shown in 
Fig. 12. 

The model is described in the following 
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a) L Written bit 


Y : Sb : InseSbso 


Fig. 13—Model of the process as a written bit cools. 


subsections. 

4.1.1 Writing 

While the film is irradiated at high power for 
a short duration, the temperature of the film 
increases rapidly, and when its temperature 
reaches the melting point, the film becomes 
liquid. After the irradiation of laser light, 
the temperature goes down. When it reaches 
the liquidus line in Fig. 12, InSb-compound 
crystals solidify at the outer part of the 
irradiated area as initial crystals. At this time, 
excess Sb atoms move from the outer part to 
the inner part of the irradiated area, as shown 
in Fig. 13. The central part still remains liquid, 
and the average composition goes towards the 
eutectic point along the liquidus line. 

During this cooling process, Sb atoms 
preserve high mobility and move quickly 
towards the center of the bit. So, after the 
temperature comes down to the melting point, 
Sb atoms form large crystals. 

4.1.2 Erasing 

Because the laser power used in the erasing 
process is low, the maximum temperature of 
the liquid film is lower than for the writing 
process. During cooling, InSb compound crystals 
first solidify at the outer part of the bits. The 
mobility of the Sb atoms, however, is low 
because of the low temperature of the bits’ 
center. Thus, Sb and InSb compound crystals 
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a) L Erased bit 


B :» 


Fig. 14—Model of the process as an erased bit cools. 
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solidify as small crystals uniformly distributed 
in the bit. 

This process is illustrated in Fig. 14. 

The authors have examined the atomic 
distribution in the bits by EPMA (electron probe 
micro analysis). The results were in good ac- 
cordance with this model”. 

The same model is applicable to the writing 
and erasing process on a rotating disk. A detailed 
analysis is now under way. 


4.2 Media characteristics 

The recording sensitivity of SeInSb and InSb 
disks is fairly good. Considering that currently 
available laser power is about 10 mW at the 
medium surface, maximum recordable velocity 
of the medium is 9 m/s, which closely cor- 
responds to the linear velocity of a 130-mm- 
diameter disk rotating at 1 200 rpm. This is 
almost satisfactory for application to low- 
transfer-rate (less than 1 Mbyte/s) devices. 
But it is not sufficient for a high-transfer-rate 
application of optical disks. For example, 
recent magnetic disks have attained data transfer 
rates of 3 Mbyte/s, with a linear velocity of 
more than 20 m/s. Thus, an approximately two- 
fold increase in recording sensitivity is necessary. 
This is a problem that must be solved in the 
future. Development of a high-power laser 
would be another way to resolve the problem. 
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Carrier-to-noise ratio of 45 dB is sufficient 
for digital recording. Stable write/erase cycle is 
on the order of 10°. The limitation is 
determined by the noise. The noise is thought to 
originate from grain growth of InSb-compound 
crystals located on both sides of the recorded 
tracks. Addition of other elements to the 
SeInSb alloy, or improvement of film fabrica- 
tion method, will resolve the problem. 

Long archaival life is an indispensable 
condition of optical memory media. The good 
chemical stability is a merit of SeInSb and InSb 
alloy films. Medium life is over ten years in 
normal environments, even if conservatively 
estimated. This is sufficient for practical use. 
It also proves that the mechanism of crystal- 
to-crystal phase change is stable. But this estima- 
tion was in terms of degradation of carrier-to- 
noise ratio. The authors are now studying 
stricter estimation by bit-error tests. 


5. Conclusion 
A new types of crystal phase change in 

SeInSb alloy film was found, and an erasable 
optical disk has been developed using this 
film. Investigations by transmission electron 
microscope strongly suggest that the phase 
change originates from a change in the distribu- 
tion of crystals, and that this change is caused 
by irradiated laser light. 

The medium characteristics can be sum- 
marized as follows. 

1) Carrier-to-noise ratio is 45 dB. 

2) Recording sensitivity is moderate; recordable 
medium velocity is 9 m/s using 10 mW laser 
power at the medium surface. 

3) Disk life is over ten years. 

4) Erasable cycles vary between 10? and 10°, 
at present. 

Though there remain some problems to be 
resolved, the disk has attained almost-practical 
characteristics. The SeInSb disk is very 
promising. The authors hope to contribute 
the progress of technology through the develop- 
ment of the SeInSb disk. 
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Erratum 


Experimental Modal Analysis of Vibrating System: A Method for Modal Parameter Identification 
and Its Software Package. Hiroyuki Kano, Fumio Matsumoto, and Junji Kaneko, FUJITSU Sci. Tech. J., 
23, 1, (1987). 

Erros found in Table 1 on page 62 of the paper should be corrected as follows: 


i a(i) b (i—1) i a(i) b(i—-1) 
5 —5.52677E-01 ee 5 5.52677E-01 
17 8.05400E-01 zak 17 8.05400E-02 


The authors wish to thank Professor S. Sugimoto of Osaka University for bringing these errors to their 
attention. 
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